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ABSTRACT 

The idea of survivahility introduced as a network parameter has led 
to so many investigations. Several measures of survivaMlity has been 
studied. The number of links and/or stations needed to be damaged to 
disrupt the system is the survivability criterion adapted in this study. 

The development of analysis procedures for directed, undirected, or 
mixed networks based on the above criterion knd use of the concepts in 
network flow and graph theory are treated in detail including computer 
program implementation of the algorithms. Finally a practical design 
algorithm for minimum-cost survivable network with respect to branch dis- 
connection using a heuristic approach and analysis techniques is described. 
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I. INTRODUCTION 



Network complexes such as railroads, highways, communication circuits, 
telephone systems, and many others which are interconnections of links and 
stations, and where there exist commodity flow, have been modelled by lin- 
ear graphs. The vertices of the graph symbolize stations or terminals 
while the branches denote lines, trunks, or circuits. The graphical struc- 
ture of such systems allows the use of topological concepts in the develop- 
ment of the theory of network. 

In 1961 the idea of survivability as a graph parameter was introduced. 
This concept vias later adapted as a network parametex'. Survivability of 
a network is defined as its capability to function as a system after dam- 
age induced by enemy attack. Survivability and vulnerability are used 
interchangeably by many authors. The tenn network reliability is more ap- 
propriately used when referring to damage caused by natural disturbances. 

The study of network survivability employs topological methods and the 
concept of network flow. There are two distinct areas, namely; determin- 
istic survivability and probabilistic survivability, A network is con- 
sidered probabilistic when the parameters involved like commodity flow and 
weights assigned to links and stations are statistical. When there is 
certainty in the existence of these parameters, the network becomes de- 
terministic. 

Up to present, network survivability has been the subject of research 
by many graph and network theoreticians. Recent works were focused on the 
mathematical formulation of survivability criteria, the development of 
analysis methods, and the synthesis of optimally invulnerable networks. 
Numerous results have already been published. 
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This thesis is based on the paper by Frank and Frisch []l0] which 
summarizes the most significant results of the researches in this area. 

The treatment of deterministic networks is restated in more detail and ex- 
tended by adapting other results in network flow theory with some modifi- 
cations to provide additional analysis procedures. To provide working 
tools for analysis and synthesisi the existing algorithms and those that 
were modified were coded for the computers. Further simplifications were 
introduced when encoding these algorithms into the programs which appear 
in the appendices attached herewith. 

In Chapter II of this report, a summary of the important elementary 
principles of network theory essential in the study of network surviva- 
bility is presented. Chapter III covers the analysis procedures based on 
some survivability criteria. The principles from which they were derived 
are also discussed. Investigation of directed, undirected, and mixed net- 
works with respect to branch disconnection and vertex disconnection are 
separately treated. 

These analysis tools may be employed in the design of optimally sur- 
vivable networks. The object of Chapter IV is design optimization by an 
heuristic approach using analysis techniques. 
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II . NETWORK FUNDAMENTALS 



2.1 NETWORK TOPOLOGY 

One of the many applications of linear graphs to engineering problems 
is the topological representation of networks. Networks may be categorized 
as "electrical networks" i "switching networks", and "flow networks". The 
latter which is often called "communication net" has associated with its 
branches some information or commodity flow. Symbols, terminology, and a 
minimum of the fundamental concepts necessary in the study of this type of 
network is discussed in this chapter. 

2.1.1 Directed, Undirected, and Mixed Networks 

A directed network is represented by a directed linear graph G=(V,r), 
consisting of a set of elements called vertices denoted by the symbol V 
and a set of ordered pairs of vertices called branches denoted by P . The 
vertex v. is called the initial vertex and the vertex v . the terminal ver- 
tex; vertices v^^ and v^ are said to be adjacent. 

A directed branch is indicated by the ordered pair (i, j). Each 
branch carries an orientation and the information flow is only in the spec- 
ified direction. 

An undirected network is represented by an undirected graph G-(V,P), 
consisting of a set of vertices V, and a set of unordered pairs of these 
vertices denoted by P . An undirected branch is indicated by []i,j3 
The flow of information is bidirectional. 

A mixed network is symbolized by a graph that has both directed and 
undirected branches. A non-oriented branch may be represented graphically 
by two parallel branches oppositely directed. 
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The network graph may have more than one branch between the same 
pair of vertices. 

2.1.2 Paths In a Network Graph 

A directed s-t path is a sequence of distinct vertices and 
branches connecting two arbitrary vertices, v and v, . It has the form 

S u 

Vs» (s,i^), V^ , (i^,!^), v^ , ... , v^ , (ij^,t), v^ 

An undirected s-t path rr^ is a sequence of distinct vertices and 
branches of the form 






[s,i^],Vi , [1^,12]. V , ... , V , [i^,t], v^ 

1 2 k 



Two directed or undirected paths tt. and rr . are branch disjoint if 

1 J 

they have no branches in common and vertex disjoint if the only vertices 

they have in common are v and v, • 

s x> 

In a sequence of vertices and branches describing a directed path, 
if (i»o) is in the path, then (i,j) is a forward branch with respect to 
the path. A branch (j,i) in the path is a backward branch with respect 
to the path. 

The network graph G is said to be connected if for any pair of ver- 
tices V. and V. there is a path n. . in G. Otherwise it is disconnected. 

1 j i-j 

The resulting subgraphs are components . 

2 . 1.3 Cutsets and Cuts 

An undirected branch cutset of an undirected, connected graph is 
the minimal set of branches, the removal of which yields a graph of two 
or more components. 

A directed branch cutset is a minimal set of branches of a directed 
graph, the removal of which disconnects all di.rected paths from at least 
one vertex of G to at least another vertex of G. 
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An undirected vertex cutset of an undirected, connected graph is 



the minimal set of vertices whose removal from G disconnects the graph. 

A directed vertex cutset is a set of vertices of a directed graph 
whose removal from G destroys all directed paths from at least one of the 
remaining vertex of G to at least another remaining vertex. 

Cutsets with respect to a vertex pair v^ and v^ are referred to as 
s-t cutsets . 

A cut (X,X) separating v and v, is a set of branches (i,j) such 

s 

that € X and v^ € X. In Fig. 2.1, if X = ^ l 






j , the set of branches (X,X) = 1(1,4), (l,3). (2,3)| 



3’ 

is a directed 



s-t cut. Branch (2,4) is oppositely directed with respect to the vertices 
s and t therefore it does not belong to the set. 




Fig, 2.1 
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2.1.4 Netv/ork Tree 



A spanning tree T = (V,U) of a connected network graph G = (V, T ) 
is a connected subgraph of G containing all of the vertices of G but no 
closed path. 

If n is the number of vertices in G, a tree contain n-1 branches. 



only one i-j path. Therefore it becomes disconnected when one branch is 
removed. 

2.2 FLOWS IN DETERMINISTIC NETWORK 
2.2.1 Branch Flow and Branch Capacity 

Aside from the topological structure of a communication net, there 
is assigned to each branch a finite weight f(i,j) or fCi,j]] called the 
branch flow from v. to v.. Also a weight c(i,j) or c[i,j] is assigned to 
each branch which is the maximum amount of commodity that can reach v. 

from v^ per unit time, called the branch capacity . 



a directed case, equations (2.2. l.l) and (2.2. 1.2) are. not necessarily 
true. 

An undirected branch can be represented by two oppositely directed 




For an undirected network. 



= -c[j,i], 
= -f[j.i]» 

f[s,t] = E f[i,j]. 



( 2 . 2 . 1 . 1 ) 
( 2 . 2 . 1 . 2 ) 
(2.2. 1.3) 



where v and v. are any chosen source and sink vertices respectively. For 
s x>. 



branches between and with capacities 
c(i.j) = c(j,i) = c[i,j]. 



(2.2. 1.4) 
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2.2.2 Flow Pattern 



The set of flows associated with the branches in G is called a flow 
pattern denoted by F, A flow pattern is said to be feasible if it satis- 
fies Equations (2.2.2. l) and (2. 2. 2. 2) for some nonnegative constant f . 

S f 1/ 

which is called the value of F. 

For all V. € V of a directed network, 

1 

r 

f , if i=s 
s,t 

0 if i/s,t (2.2.2. 1) 

-f . if i=t 

^ s,t 

0 < f(i,j) < c(i,j) for all i and j (2. 2. 2. 2) 

This means that the net flow out of the source v is f . and the 

s s, t 

flow out of the sink v, is -f . i whereas the net flow out of an intermed- 

s I 

iate vertex is zero. This latter relation is also true for an undirected 
network. 

A branch (i,j) is saturated if f(i,j) = c(i,j) and unsaturated 
otherwise. It is said to be empty if f(i,j) = 0. 

2.2.3 Residual Capacity and Augmentation Path 

The residual capacity i'(n , ) is defined both for directed and un- 

Sit/ 

directed network as 

r(n ,) = Min Min rf(i,j)1 , Min f c(a,b) - f(a,b) 

' L backward forward 

branches brcinches 

(i.j) € (a,b) € (2.2.3. l) 

A path IT , is saturated if r(n . ) = 0 and unsaturated otherwise. 

S I S 1 1/ 

Additional flow can be sent through an unsaturated path obtaining a new 
flow pattern. To achieve this, add r(n . ) to all forward flows and 

S I u 

subtract r(n .) from all backward flows along the path. 

S I 

If a path between two vertices has 
f(i,o) < c(i,j) 




f(i,V) - f(V,i) = 
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aaid 



f(j.i) > 0, 



then we call this a flow augmentation path . The flow along this path cam 
be increased according to the previously stated rule. 

2.2.^ Maximizing Network Flow 

One of the most important and basic theorems in network flows is 
the "Max- flow Min-cut theorem" formulated by Ford and Fulkerson 1^7] . The 
problem of obtaining the maximum obtainable flow value f , through a net 

S I X 

with reference to two vertices v ajid v, is solved using this theorem. 

S X 

Max- flow Min-cut Theorem: The maximal flow value obtainable in a 

network G form v to v, is the minimum of the value of a cutset taken over 
s t 

all cutsets separating v^ and v^. 

To state the theorem mathematically, let: 

5* , be the value of the minimum s-t cut 

S j X 

c(X^,X^) be the capacity of a cut (X^,X^) 

•j (X. ,X. ) I be the set of all s-t cuts in the net 

1 1 s , t. J 



then: 



Max [f > Min [c(X ,X ) ] = . {2.2.k.l) 

The value v of the cut (X,X) is defined as 

V = f(X,X) - f(X,X) < c(X,X) {2.2.h.2) 

The Ford and Fulkerson Labeling and Augmentation algorithm [[ 7 ] based 
on the above theorem provides a method of increasing f , systematically 

S J X 

from an arbitrary feasible flow pattern with value less than ^ , and ter- 

S J X 

rainates when f , - T +• resulting flow pattern is maximal. 

S J X S J X 

The algorithm will be stated in a modified form in Chapter III. 
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2.3 MATRICES ASSOCIATED WITH SUEVIVABLE NETWORKS 

Network graphs may he structurally drawn or he descrihed. hy matrices. 
The latter representation is useful for computer calculations. This sec- 
tion enumerates several matrices used in the study of survivahle networks. 
2.3tl The Incidence Matrix 

Let G = (V, r* ) he the directed graph of the network consisting of 

n vertices and m branches. Let A = [a. he an n x m matrix whose i-jth 

entry is a. .. Then A is the incidence matrix of G if 
1 > J 

+1 if a branch is directed away from v^ 



a. . = < 

ii J 



-1 if a branch is directed towards v. 

1 

0 otherwise 



If G is an undirected graph, A = []a. .]] is the incidence matrix if 

r 



a. . = < 

i.J 



1 if branch is incident at v. 

1 



0 otherwise 



2 . 3.2 Connection Matrix 

An alternate representation of a non-oriented network is the con- 
nection matrix K = [|k. K is an n x n symmetric matrix whose i-jth 

1 > J 



entry is k. . where 
1 1 J 






k. . 

ii J 



1 if [i.j] € P 






0 otherwise 



For a directed network, a -1 is entered for an oppositely directed branch. 
2 . 3.3 Branch- capacity and Terminal Capacity Matrices 

The other way of representing a communication net of n vertices is 
by using a symmetric matrix C = [c. .J called the branch- capacity matrix 

1 1 J 

where c. . = c. . is the sum of the capacities of the branches connected 

J.i 

between v. and v.. For an oriented network, the C-matrix is asymmetrical. 
^ J 
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For any two vertices in G, there exist a finite, maximum obtainable 



capacity of transferring information from one vertex to the other called 

the terminal capacity. The terminal capacity matrix T = [t. represents 

1 » J 

the maximum flow that can exist between all the possible vertex pairs of 

G, The T-matrlx is also asymmetrical for a directed network, 

2.3*^ Connectivity and Redundancy Matrices 

The connectivity matrix is essentially a T-matrix with each branch 

of G or G assigned a maximum capacity of unity. 

An n X n matrix R = []r. .]] is a redundancy matrix where the i-jth 

1 » 0 

entry r^ ^ is the number of branch/vertex disjoint paths between v^ and 
Vy It is symmetrical for a non-orlented net. A uniform R-matrix is one 
which has Identical entries except the main diagonal. 

The main diagonals of these matrices have infinity as entries how- 
ever for computer calculation zero entries are used. 



2.4 DISCONNECTING SETS OF A NETWORK 

A connected network represented by G can be disconnected by deleting 
certain branches or vertices. 

An s-t branch cutset is a set of branches whose deletion breaks all 

directed paths from v to v , ; likewise an s-t vertex cutset is a set of 

s t 

vertices' whose deletion breaks all directed paths from v^ to v^. A mixed 
s-t cutset is composed of branches and vertices. 

Clearly if the capacity of each branch or vertex is set to unity, 
then the value of the cutset is numerically equal to the number of com- 
ponents of the cutset. 

Let a .be the number of elements in the smallest s-t vertex cut- 
s, t 

set; 

, be the number of elements in the smallest s-t branch cutset; 

S f ’b 
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be the number of elements in the smallest s-t mixed cutset; 



then the disconnecting sets of the connected graph aj:e as follows: 

W = Min [ ,] (2.4. 1.1) 

S j X- 

s,t € V 
Sj?^t 

T^MinCr .] (2. 4. 1.2) 

S j X 

s,t € V 

s/t 

<r=Min[0'g^J _ (2. 4.1. 3) 

s,t € V 

s^t 



The smallest disconnecting set of a connected graph denoted by v is 
defined as 



V = Min [ (d , r , <r ] 

However it was shown in [[3]] that 
V Min [ hi 

S > X 

s,t 6 V 

s/t 



(2. 4.1. 4) 

(2. 4.1. 5) 
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III. ANALYSIS OF DETERMINISTIC SURVIVABLE NETWORKS 



3.1 SURVIVABILITY CRITERIA 

A significant criterion in network analysis is its survivability or 
vulnerability. In a coimunl cation net, this parameter refers to the de- 
gree in which the system remain functional after an enemy attack or nat- 
ural disaster. 

No general survlvable criteria has yet been formulated however several 
measures of survivability has been investigated namely s 

1. The number of links and/or stations that must be destroyed before 
communication is disrupted 

2. the minimum number of links that must be deleted in order to iso- 
late a group of stations from communication with the remainder of the net 

w, 

3. the minimum number of isolated subnets that would result from an 
attack or destruction. This is termed the "independence number" of the 
network [l], 

4. fraction of stations which survive an attack [ 2 ], [6^, 

5 . fraction of stations that can be reached from a given point by 
a path of no more than r links after an attack j^8]. 

The first three criteria apply to deterministic network while the 
last two criteria are associated with probabilistic net. 

Analysis of large networks based on the first criterion is treated 
in this thesis. This index of survivability suggests the application of 
the concept of disconnecting sets. More specifically, the parameter U 
is a measure of vulnerability to station destruction and is a measure 
to link destruction. 
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3.2 ANALYSIS WITH RESPECT TO BRANCH DISCONNECTION 
3.2.1 s-t Branch Connectivity 

The problem discussed here is the finding of the maximum number 

of branch disjoint s-t paths between any given pair of vertices v and 

s 

v^ in G. Setting the branch capacities to unity and applying the max- 

flow min-cut theorem, the maximal flow f . from v to v, is obtained 

s,t s t 

which is numerically equal to the smallest brsmch cutset . . This 

S > X- 

value further represents the maximum branch disjoint paths between the 
given vertex-pair. 

The Ford-Fulkerson algoirithm solves a general max-flow problem be- 
tween two vertices, the source and sink. It is stated here in simplified 
form to adapt it to the special case being investigated. 

A few terms are defined below before stating the algorithm: 

1. A vertex is labeled by an ordered pair (i,+). For a particu- 
lar vertex Vy the index i denotes the preceeding labeled vertex v^. 

The second teim in the label is plus if for a forward branch, f(i,j) = 0 
and a minus sign if for a backward branch, f(j,i) > 0. Otherwise a ver- 
tex is unlabeled. 

2. A particular vertex v^ is scanned if its adjacent vertices v^ 
are labeled from it. Otherwise v is unscanned, 

3. A vertex is considered to be in one of the three states, namely 

a. unlabeled and unscanned 

b. labeled and unscanned 

c. labeled and scanned 

The follo:n.ng steps determines the s-t branch connectivity of a 
directed or mixed network: 

1, Represent the network by a directed linear graph. Number the 
vertices consecutively from 1 to n. 
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2. For each branch (i,j) € P , assign c(i,j) = 1 and = 0. 

3. Choose the source and sink vertices v and v. . 

S w 

4. Apply Algorithm 1 

Algorithm 1 (Modified Ford and Ptilkerson Algorithm) 

Labeling Routine ; 

a. Erase all vertex labels 

b. Label v by (0,+), Source is now labeled and unscanned, 

s 

all other vertices are unlabeled and unscanned, 

c. Select any labeled unscanned vertex v. (including v ). 

1 s 

Assign labels to all unlabeled vertices v^ adjacent to v^ such that 

l) for forward branch (i,j) with flow f(i,j) = 0, label 

V. by (i,+). Vertex v. is now labeled and unscanned. 

Z) for backward branch (j,i) with flow f(j,i) > 0, label 

V. by (i,-). Vertex v. is now labeled and unscanned, 

J 3 

With all adjacent vertices to v^ labeled and unscanned, 

V. is now scanned, 

1 

d. Repeat Step c until either ' . 

l) is labeled or 

Z) no more labels can be assigned and is unlabeled. 

In Step d(l), an augmentation path is found and an increase 
of flow is possible. Proceed to Augmentation Routine. 

In Step d(2), the algorithm terminates. Flow is now max- 
imized. The disconnecting cutset are the branches joining the labeled 
and unlabeled vertices. The value of the cut is the branch connectivity. 
Augmentation Routine ; 

a. Let z = t and go to Step b, 

b. If label on v is (q,+), set f(q,z) =1, If label on v 

z z 

is (q,-), set f(z,q) = 0. 
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c. If q = s return to Step a of the Labeling Routine. Otherwise 
set z = q and go to Step b of the Augmentation Routine. 

Example ; A mixed network 

Fig. 3.1(a) is the linear graph representation of the network. 

In (b), the non-oriented branch (3>6) is split into two oppositely direct- 
ed branches (3i6) and (6,3). The first labeling obtained the augmentation 
path s, (s,3) i3i (3i?) i7i (7it) . In (c) the second labeling found the aug- 
mentation path s, (s,4) ,4, ( 4 , 3 ) ,3i (3i6) ,6, (6, t) . In (d) the sink vertex 
cannot be labeled since (6,8) and ( 7 , 8 ) are already saturated and (5i8) 
is a backward branch with zero flow. The algorithm terminates. The cut 
is ^ (6,8) ,(7i8) I and s-t connectivity is two. 

Flov; chart for this analysis is in Appendix A-1 and the computer 
program MAXCON-IA is in Appendix B-2. To minimize computer storage, a 
vertex label is further modi.fi ed by using a single signed integer (+i) 
instead of the ordered pair (i,+). V/hen i is negative the magnitude is 
obtained by just multiplying by minus 1. 

In the program, vertex scanning is accomplished by using two stor- 
age vectors. The first vector is a list of the labeled vertices being 
scanned while the second provides storage for the vertices being labeled. 
Once the scanning process of all the vertices listed in the first vector 
is completed, the contents of the second is transferred to the first which 
in turn becomes the list of vertices to be scanned. 

Algorithm 1 is applied to the analysis of undirected network with 
slight modification. For a non-oriented net, each branch flow is con- 
sidered bidirectional with a maximum capacity of unity. However, in no 
case can a positive back-flow occur since each augmentation path has as 
its elements all forward directed branches with respect to v and v. . 
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Fig. 3.1 
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Consequently there csinnot be any negatively labeled vertex. Therefore the 
tests for a backward branch and a negative-labeled vertex, and the adjust- 
ment of backward flow are not necessary in this case. 

These steps are deleted in Algorithm 1 and MAXC0N-2A (Appendix 
B-4) is developed for determining s-t branch connectivity of an undirected 
network. 

3 . 2.2 Multiterminal Branch Connectivity 

For an n-statlon network, its graph would have n(n-l)/2 possible 
vertex pairs. To apply Algorithm 1 n(n-l)/2 times in order to obtain the 
connectivity between all the possible vertex-pairs would be lengthy. 

Gomory and Hu have shown that for an undirected network with 

n vertices, any flow in the network is numerically equal to some flow in 
a maximal spanning tree. Since there are only n-1 branches in a spanning 
tree, then there can only be n-1 numerically different flows possible. 
Therefore all the n(n-l)/2 possible maximal flows that can exist in the 
network can be deduced after doing n-1 flow maximizations. 

Applying the above principle and the max-flow min-cut theorem in 
the analysis of non-oriented nets, the multi-terminal network connectivity 
problem can be solved 'by only n-1 separate applications of Algorithm 1. 

The analysis procedure is embodied in the following algorithm. 

Algorithm 2 

1. Choose two vertices v and v. in the undirected graph G. 

S Xf 

With branch capacities set to unity, find an initial s-t branch cut using 
Algorithm 1. Represent this cut by a generalized vertex tree with con- 
densed vertices X and X connected by a branch with capacity equal to the 
computed cut value T .. 

S J tr 

2. Draw a condensed network graph G^ where X is represented 
by a single vertex. Select two vertices Vj^ and v^ in X and find the 
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smallest k-1 cut by Algorithm 1. Form the generalized vertex tree by 
dividing X into and X^ vdth ^ as the capacity of the branch con- 
necting them. In this tree, X is adjacent to X^ if it is in the same 
side of the cut as X^, or X is adjacent to X^ if it is in the same side 

of the cut as X^, 

3 . The above procedure is repeated. With each step, a tree 
of generalized verrbices connected by branches vdth corresponding cut 
values is obtained. To proceed vdth the computation, 

a. Select a generalized vertex X^ 

b. Condense each component in G except X^ into a single 
vertex forming the condensed netvjork graph G^, 

c. Choose two original vertices v and v, in X. as source 

S t 1 

and sink and apply Algorithm 1 to G. to get 'J* . , 

X S > X 

The cut obtained in this process splits X^ into two parts 

X. and X. , This is represented in the generalized vertex tree by re- 
^1 ^2 

placing X. by the condensed vertices X. and X. connected by a branch 
bearing the cut value. All other branches and vertices in the tree are 
unchanged except those which were formerly connected to X^. Apply the 
adjacency mle as in Step 2. 

4. Repeat the process until these generalized vertices con- 
sist of exactly one vertex each. This results in a spanning tree. Thus 
for any chosen v and v. , the s-t connectivity is the smallest branch 
weight in the unique s-t path in the tree. 

An example for the application of the above procedure is given 
in Fig. 3«2. The generalized vertex tree and the condensed network G^ 
resulting in each flow calculation are indicated. For the 12-vertex net- 
work, the spanning tree is formed after 11 max-flow computations. 



23 



Fig. 3.2 
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A flowchart for this analysis procedure is given in Appendix A-2, 
This is further encoded into a computer program entitled HAXCON-2 appear- 
ing in Appendix B-3. In this program the printout is not a spanning tree 
but rather a symmetrical branch-connectivity matrix which is essentially 
a terminal capacity matrix. Its entries are derived from the cut values 
appearing as branch weights of the trees generated after every max-flow 
calculation where a cut isolates a single vertex from the rest of the con- 
densed network graph G^. The connectivity between this single generalized 
vertex and all the rest of the vertices in G are the updating entries. 

Other information relevant to HAXGON-2 are stated in the program 
docum ent at i on , 

The concept of spanning tree cannot be applied to the analysis of 
a directed or mixed network as generally these networks are not symmet- 
trical. To obtain a multiterminal branch-connectivity matrix for an ori- 
ented net, it is necessary to do n(n-l)/2 flow maximizations by Algorithm 
1. Up to present, there is no available method that could reduce this 
number of calculations however Frisch has developed the Flow Variation 
Algorithm C 12 ] which effectively reduces the number of steps required to 
maximize f 

s, t 

■ The ideas underlying this new algorithm are as follows: let v , 

s 

V , , and V. be three arbitrary vertices in a graph G, Rather than apply 

S X> 

Algorithm 1 to independently maximize f and f . , first maximize f , . 
by Algorithm 1 then maximize f^ ^ in the usual manner until a vertex v^ 
is labeled for which f(g,V) > 0. Then use the flow pattern D(s' ,t) to 
find F(s,t) , 

Applying this basic idea to solve the multiterminal connectivity 
problem which is a special case of a general flow problem, all that is 
necessary is to slightly modify the algorithm to fit into the problem. 
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The result is a simplified version of the original algorithm. The following 
notations are used; 



1 

-p 


flow path for vertex-pair v and v. 

S If 


^s' |t 


flow path for vertex-pair v , and v. 

S Xr 


F(s,t) - 


‘ flow pattern for vertex-pair v and v. 

S X 


D(s’,t)- 


flow pattern for vertex-pair v , and v. 

S X 


f(i.j) - 


■ a branch flow in path n . 

S 1 X 


<i(i.j) - 


a branch flow in path n , , 

S 1 X 


A 

TT 

S|g 


section of path rr , from v to v called "Augmenta- 
S 1 1 . s g 

tion Path" 


T 
TT 1 

s' ,g 


section of path tt , . from v , to v called "Trunca- 
s 1 r s g 

tion Path" 


H 


• section of path n , , from v to v, called "Homing 

s , t. g t 

Path" 


- 

m,m 


• a loop detected by the Homing Routine; m is the first 
vertex in the loop 


- 

P.P 


a loop detected by the Truncation Routine; p is the 
first vertex in the loop 



The subscript c whenever it appears with the symbols indicates the 
th 

c stage of max- flow calculation. 

The following routines are employed in the algorithm: 

Routine A - Detecting an Augmentation Path 

1. Label v by (s,+). v is now labeled and unscanned. All 

s s 

other vertices are unlabeled and unscanned, 

2, Select any labeled and unscanned vertex v^. Assign labels 
to all unlabeled vertices adjacent to v^ such that 
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a. for forward branch (i, j) € T with flow f(i,j) = 0, 

label V. by (i,+). v. is now labeled and unscanned, 

J 0 

b. for backward branch (j,i) € P with flow f(j,i) > 0, 
label V. by (i,-), v . is now labeled and unscanned. 

With all adjacent vertices to labeled and unscanned, 

is now scanned. 



3, Repeat Step 2 until either: 

a, labeled or 

b, a vertex v is labeled with d (g,V) > 0 

c, no more vertices can be labeled and v^ is unlabeled. 
Observe that except for Step 3(l3) this routine is the Ford 

and Fulkerson Labeling Routine, 

Routine T - Detecting a Truncation Path 

1, Let z = g 

2, Find one vertex v. such that d (i,z) >0 and label v. by 

1 c 1 

(z) 

3, Let z = i 

a. If V has assumed a value i before, terminate, 

z 

b. If i = s' , terminate, 

c. Otherwise return to Step 2, 

T 

If the routine terminates in Step 3(a) » a loop L has been 

p,p 

detected rather than a path. 

Routine H - Detecting a Homing Path ■ 

1, Let q = g 

2, Find one vertex v. such that d(q,j) > 0 and label v. by 

(q) 



3, Let q = j 

a. If q has assumed the value j before, terminate. 
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'«||0 



I 



b. If j = t, terminate. 



c. Otherwise return to Step 2. 

H 

If the routine terminates in Step 3(a), a loop L has been 

' m,m 

detected instead of a path. 

The main algorithm for obtaining the multiterminal branch-connec- 
tivity of an oriented network is now stated in terms of Routines A, T, 
and H. 



Algorithm 3 



1. Choose two vertices v and v. . Using Routine A, find an 

S "C 



initial s-t cut. 



2. Using Routine A find an s-g augmentation path tt or an 

s , g 

s-t augmentation path n . . In the former case go to Step 3 . In the lat- 

S I Xr 

ter case go to Step 8, If no such path can be found, terminate. 



H 



3 . Using Routine H find a Homing Path n . or a Homing Loop 

L . In the former case go to Step 4; in the latter case go to Step 
m,m 

5(a). 

T 

4. Using Routine T find a Truncation Path rr , , or a Trunca- 

S’ I x» 

T 

tion Loop L . In the former case go to Step 6; in the latter case go 
P>P 

to Step 5(b) . 



5 . a. Set to zero all d^(i,j) for all branches (i,j) in the 



.H 



loop L . Erase all vertex labels and go to 2. 



m,m 



b. Set to zero all d (i,j) for all branches (i,j) in the 



loop L . Erase all vertex labels and go to 2. 
P.P 



H 



6, Trace and erase g-t path. For every branch (i,j) in n , , 



set 



dc+i(i,j) = 0 
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7. 



T 

Truncate the s'-g path. For every branch (i,j) in n , , 

s »g 

set = 0. 

8. Augment the s-g path. For every branch (i,j) in rr'^ , set 

Sf g 

f^+l(l.j) - 1. 

Erase vertex labels and return to Step 2. 

If a "loop of flow" is detected, the flow in this loop must be 

set to zero as in Step 5 a-s otherwise a nonfeasible flow will be achieved. 

As an example, consider the graph in Fig, 3* 3* bet the Indicated 

flow pattern D(s' ,t) be the result of initial maximization of f , , . 

Designate v^ as the source for the second stage of maximization. Step 2 

of the algorithm labels v^ and detects it as a vertex v^. Step 3 finds 

the sequence 3»(3»6) ,6,(6,8) ,8 as the Homing Path n .. The path s',(l,3)> 

St'^ 

T A 

3 is the Truncation Path n , , and s,(2,3)»3 the Augmentation Path n 

s > g s »g 

step 6 sets to zero the hranch flows d (i|j) from to v, and replaces 

each branch flow in the path by b^^(i,j) - !• Step 7 truncates s’ -3 path 

by setting d^(l,3) to zero, and Step 8 augments s-3 path by setting 

f^_j^^(2,3) ■■ !• One s-t path is now processed and the algorithm will 

search for csnother one until no more augmentation paths can be found. By 

then, f , will be maximal, 
s, t 

Flowchart of this algorithm is in Appendix A-3 and its computer 
implementation MAXCON-1 is presented in Appendix B-1, 
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3.3 ANALYSIS WITH RESPECT TO VERTEX DISCONNECTION 



The criterion of survivability attached to this problem is the number 
of stations that must be deleted from the network in order to disrupt com- 
munication. This is the smallest disconnecting set W, of a network 

graph. Given any pair of vertices v and v . , assign unit capacities to 

s 

all vertices except v and v. (which have infinite capacities), then the 

S 

smallest s-t vertex cutset W + is numerically equal to the maximiim num- 

S f 

ber of vertex disjoint s-t paths. The value of («J .is the network ver- 

S f 

tex-connectivity with respect to the given source and sink. 

3.3*1 s-t Vertex Connectivity 

To solve the problem of finding W . , each vertex in G is split 

S f *0 

into tvjo vertices connected by a branch of unit capacity forming a new 
graph G'. Algorithm 1 is then applied to G'. The result is a branch 
cutset composed only of the unit capacity branches. The regular network 
branches are considered having infinite capacities. This approach is 
however very inefficient since the number of vertices are almost doubled. 

An algorithm known as .Vertex-Pair Connectivity derived from Algo- 
rithm 1 was developed by Frisch (^ll^ which applies directly to graphs of 
directed, undirected, or mixed networks without splitting the vertices. 

The vertex-splitting process is absorbed in the modification of the label- 
ing routine of Algorithm 1. The algorithm is stated here in a different 
form. Several terms and rules must first be defined: 

1. A vertex is labeled with an ordered triple (i,+»M). The 
first and second entries are defined as in the original Labeling Routine. 
The third entry is a zero or one as will be explained in the algorithm. 

2. A vertex is strongly labeled if the second entry is a minus, 
and weaJcly labeled if it is a plus, 

3. A circle around the second entry denotes the vertex is 



seamed. 



32 







4 



’V 





If a vertex is weakly labeled and unscanned, and then 
receives a strong label, the weak label Is erased. However if the label 
is scanned, it should not be erased. 

Algorithm 4 
Labeling Routine 

1. Erase vertex labels. Label v by (s,+,0). v is now 

s s 

labeled and unscanned. All other vertices are unlabeled and unscanned, 

2, Select any labeled unscanned vertex v^. If none, ter- 
minate. 



a. Suppose v^ is weakly labeled: 

1) If f(i,j) = 0, (i,j) € P and j / s, weakly 

label all vertex v. by (i,+,0). v. is now weakly labeled and unscanned. 

0 1 

2) If f(j,i) = 1, (j,i) € r , j / s, and v. not 

strongly labeled, strongly label v. by (i,-,0). v. is now strongly label- 

0 0 

ed and unscanned. 



Scan v^ by encircling the + in the v^ label, 
b. Suppose v^ is strongly labeled 

1) If f(i,j) = 0, (i,j) € r , and j ^ s, weakly 



label all unlabeled vertices v. by (i,+,M). v. is now weakly labeled and 

0 0 

unscanned. 



2) If f(o,i) = 1, (j,i) € r , i / s, and v. not 

strongly labeled, strongly label v. by (i,-,M). If v. is weakly labeled, 

0 J 

erase the weak label, v. is now strongly labeled and unscanned. 

0 

In (b), if v^ is both weakly and strongly labeled, 
set M = 1, otherwise M = 0, Scan v^ by encircling the minus sign in the 
v^ label, 

3. If v^ is labeled, proceed to Augmentation Routine, 
Otherwise return to Step 2. 
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Au^entation Routine ; 



1. Let z = t 

2, a. If is weakly labeled by (q,+,M), set f(q,z) = 1. 

b. If V is strongly labeled by (q,-,M), set f(z,q) = 0. 
z 

3. If q = s return to Step 1 of Labeling Routine. Otherwise 
proceed to Step 4, 

4, Let z = q. If M = 1 return to Step 2(b), Otherwise return 
to Step 2(a) o 

To identify the s-t vertex cutset, let U be the set of all unlabeled 

vertices, and L be the set of all labeled vertices. Further let L be the 

u 

set of labeled vertices connected to an unlabeled vertex by a branch di- 
rected from a labeled vertex to an unlabeled one, then the s-t vertex cut 
is given by; 

= L^ if s ^ L^ (3.3. 1.1) 

\ ^ - ^S^ ^ ^ ^ (3.3. 1.2) 

Since the source vertex cannot be a component of the vertex cut, Eq. 

( 3 . 3. 1.2) shows that if v^ is connected to unlabeled vertices v^^, then 
the vertices v^^ are some component of the cut. 

The graph in Fig. 3 .^ illustrates the application of the algorithm 
for a directed network. Assume that the path indicated by solid arrows 
is a result of initial labeling and augmentation. Assume further that 
the second application of labeling routine is just completed and v^ is 
labeled as indicated by the vertex labels. When the augmentation routine 
is applied, it will backtrack the path in the usual manner, however, on 
reaching v^ it has to select the strong label since v^ has M = 1 . If the 
weak label were selected, a non-feasible flow would result. This explains 
the purpose of the third entry in the label. The only case where two 
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scanned labels can occur on a vertex is the situation in where a flow 
of unit magnitude already exists through the vertex and this vertex re- 
ceives labels in the succeeding iteration of the algorithm. If the weak 
label in v^^ was unscanned when it received a strong label then the weak 
label should have been erased and v^ would have been labeled (4,+,0). 

Program MAXC0N-3A in Appendix B-6 is a computer implementation of 
Algorithm 4. A different way of circumventing the situation in vertex 
of Fig, 3*^ was applied. For a vertex like v^^^ where a flow already 
exists, if in the succeeding labeling process v^,^ receives a weak label, 
it will not be scanned, but if it receives a strong label it will be scan- 
ned. The latter vjould indicate that the existing outward flow from v^^ can 
in effect be diverted to another path. In the former, it would be waste- 
ful to continue seeking another path passing through v^^ when the vertex 
is already saturated. In no case then can a vertex acquire two scanned 
labels. The third entry in the label is therefore no longer necessary. 

The same form of label as in MAXCON-IA is used in this program. 

All steps in the algorithm that are associated with the third entry in 
the label are deleted, 

A flowchart of Algorithm 4 is in Appendix A-4, 

3.3*2 Multiterminal Vertex Connectivity 

The Flow Variation Algorithm of Frisch as presented in Sec. 3*2 
maximizes flows in a finitely branch-weighted network, and it was readily 
adapted to solving the multiterminal branch-connectivity. In this sec- 
tion the problem is different in the sense that the vertices control the 
information flow. 

In order to utilize the principle of the Flow Variation Algorithm 
here, the Routine A of the algorithm must be altered. Recall that Routine 
A is the Ford-Fulkerson Labeling Routine. If this routine is replaced by 
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the Labeling Routine of Algorithm 4 which is Frisch's Vertex-Pair 
Connectivity Algorithm, it would enable the Flow Variation Algorithm to 
handle the multiterminal vertex-connectivity problem. Consequently the 
method in which the s-g and s-t paths are augmented must also conform 
with the Augmentation Routine of Algorithm 4. 

With these changes, the algorithm flowchart is the same as in Ap- 
pendix A-3, however there is a significant change in the computer pro- 
gram MAXCON-3 as presented in Appendix B-5 compared to MAXCON-1. 

The solution given earlier for obtaining the multi terminal vertex- 
connectivity of a directed network applies equally well for undirected 
networks. Sometimes one finds it impractical to compute all the entries 
in the entire connectivity matrix especially when investigating a very 
large network. It is often sufficient to just verify whether a network 
satisfies a certain minimum degree of connectivity. 

V/ith this analysis criterion it is only necessary to verify that 
the network is at least R-connected. The connectivity being referred to 
here is the vertex-connectivity which we know is the smallest. 

The following theorems of Kleitman [l43 greatly reduce the number 
of verifications necessary rather than doing n(n-l)/2 separate tests for 
the existence of R vertex-disjoint paths. 

Theorem: 

To verify that all n(n-l)/2 vertex pairs in G can be connected 
by at least R vert ex- disjoint paths, it is only necessary to verify that 
some set of R vertices are each connected by at least R vert ex- disjoint 
paths with every vertex in the graph. 

Two other theorems resulting from the above are stated here in an 
algorithm form and are in themselves the analysis procedure. 
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Algorithm 5 



1. Choose ajiy vertex, say as source. 
Let i = 1, 



other vertices 



a. Select a vertex v. in G, as sink. Remove from G. all 

J 1 1 

k = 1,,.,,1 previously designated as sink that 



{\1 ' 



are adjacent to v.. If 1 > R go to Step l(b). 

J 

l) For V. connected to v. 

0 1 

Excluding the direct path i-j, use Algorithm 4 to 
verify existence of q vertex-disjoint paths from v^ to v^. Set s = q + 1. 

If s < R go to Step 4. Otherwise go to Step 1(b). 

b. Are all n-i vertices processed? If not go to Step 1(a), 
otherwise go to Step 2, 

2. Decrease R by i. If R = 0, terminate! otherwise go to 

Step 3* 

3. Remove v. from G. , Choose vertex v. , . as source. Increase 

1 1 l"rl 

i by one. Return to Step 1(a). 

4. Record s-connectivity between v^ and Vy Terminate. Net- 
work is not R-connected, 

The first theorem may be restated in a more specific sense; 

If there are at least R vertex-disjoint paths between v^ and 
v^ and v^ and v^ in G, then to verify the existence of R such paths be- 
tween v^ and v^ in G, we need only look for R - 1 such paths in G - • 

For the proof of these theorems and illustrative example, the read- 
er is referred to the original paper by Kleitman, 

A flowchart for computer implementation is given in Appendix A-5. 
Program MAXCON-4 in Appendix B-7 is encoded to analyze an undirected net- 
work for a uniform R-connectivity where the value of R is initially 
specified. 
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IV. DESIGN OF LOW-COST SURVIVABLE NETWORK 



4. 1 INTRODUCTION 

There are now a number of synthesis procedures for svirA/ivable networks 
appearing in literatures elsewhere. These are summarized by Frank and 
Frisch in [[lO^ and are categorized according to link or station damage. 

The synthesis approach are all aimed towards maximization of redundancy 
between stations and satisfying network redundancy requirement with the 
assumption that cost is either uniform for each link/station or not con- 
sidered at all. 

In the design of practical flow networks, the cost of construction 
increases with increase in redundancy therefore the design approach must 
take into account non-uniform cost for the construction of links or sta- 
tions, The goal is a minimum cost network that satisfies a prespecified 
redundancy requirement. So far there is no exact analytical method that 
could be applied to solve the problem. A procedure to search for a low- 
cost network with respect to vertex disconnection is presented by Steig- 
litz, etal [l8[] using an heuristic approach. The algorithms they devel- 
oped were based on the techniques for an approximate solution to the 
traveling salesman problem. The salesman's tour is a Hamiltonian circuit 
which is a special case of the flow network having a redundancy of two. 

4.2 DESIGN WITH BRANCH DISCONNECTION CRITERION 

In some cases branch destruction is more likely than vertex destruc- 
tion and it is just as necessary to obtain a design method using the 
former criterion as the latter, A procedure is presented in this chapter 
for design of a low cost network (with high probability of being optimal) 
with at least r. . branch disjoint paths. 
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The same procedure used by Steiglitz, etal is employed consisting of 
two parts. The first part called the Starting Routine generates a fea- 
sible initial network (one which satisfies redundancy requirement) and 
the second part is the Optimizing Routine which maJces an iterative im- 
provement on the network cost by branch exchanges while maintaining fea- 
sibility, Portions of the procedure that differ from that of Steiglitz 
are the feasibility testing of the network and the optimizing routine. 

The method of optimization that will be presented is an adaptation of the 
recently published heuristic algorithm for the traveling salesman problem 
by Lin and Kernighan Cl6]. 

The heuristic design algorithm is stated as follows; 

Algorithm 6 

1. Generate a feasible initial network G that satisfies the 
specified redundancy R using the Starting Routine. 

2, Attempt to find an improved feasible solution G' by some 
local transformation using the Optimizing Routine. 

3. If a similar network with a lower cost is found, then re- 
place G by G* and repeat from Step 2. 

4, If no local improvement can be found, G is a locally optimal 
solution. Repeat from Step 1 until computation time runs out or the so- 
lution is satisfactory. 

4.2.1 The Starting Routine 

This routine was developed by Steiglitz, etal in [l8^, 

1, Plot vertices representing stations in a Euclidian 2-space. 

2. Number the vertices uniformly at random, 

3» Assume cost to be proportional to distance and obtain the 
cost matrix from the Euclidian distances between vertices truncated to an 
integer. 



.. j 




4. Assign redundancy to each vertex. 

5 . Add branches to the network one at a time. Reduce by unity 
the assigned vertex requirement every time a link is connected to that 
particular vertex. A branch is added in this manners 

a. Connect the vertex with highest updated requirement to 
another vertex with the next highest updated requirement; for all vertices 
having equal requirement, choose one that results in lower cost. 

b. In cases of ties in requirement and cost, these are 
resolved by choosing the vertex highest on the list. 

c. No parallel branches are allowed. 

6. When the requirement in all the vertices are reduced to 
zero, test feasibility of the initial solution using Algorithm 2. 

7 . If solution is feasible go to Optimizing Routine, otherwise 
repeat from Step 2. 

This routine is made nondeterminlstic by numbering the vertices 
uniformly at random at each start. Randomized starting solution has been 
observed experimentally to produce a variety of starting networks as con- 
trasted to constructive solutions which are deterministic where more than 
one solution may not be obtained. 

The analysis tools presented in Chapter III come in handy for test- 
ing feasibility. For the initial network, the one that fits is Algorithm 
2 which determines the overall branch connectivity of the network. A 
complete checkout need not be carried any further once infeasibility is 
detected. 

4.2.2 Optimizing Routine 

The optimizing procedure applied by Steiglitz is the X-change where 
two branches are interchanged at a time, and determined whether the trans- 
formation is favorable and preserves feasibility. Lin and Kemighan have 
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lately experimented on a more systematic and effective way of branch 
transformation in a network with 2-redundancy. Their results are applied 
here for a network with redundancy greater than two. 



The basic ideas underlying the transformation are as follows; 
Consider a nonoptimal but feasible netvjork G with cost c(g) and an opti- 
mal or near optimal network G' with the same number of vertices and with 
cost c(G') < c(G). G is nonoptimal because it has k links that 

are "out of place". To make G optimal they "should be replaced by k links 

find k by identifying sequentially, element 

by element, two disjoint sets of links X = {'1 \ 

using a gain criterion. Let the costs of and be |x^| and |y^j 
respectively. Define = I ^i I 1 ^il gain from exchanging x^ 

with y. . For some i, g. may be negative however we want the final gain 
^ k^ 

to be positive, i.e.. Eg. > 0. We are looking for sequences of g. ' s whose 

i ^ k ^ 

partial sum is always positive. When Eg. <0, the search is stopped. 

1 

The k linlcs in X are deleted or "broken" and replaced by k links 
in Y resulting in G* which has lower cost. The process is iterated until 
no further reduction in cost is attained. 

We now state the routine more formally assuming that a starting 
solution has been generated: 

1. Set P* = 0. (P* is the best improvement so far). Choose 

any vertex as t^ and let x^ be one of the edges of G incident to t^. Let 
i = 1. 



2. Designate the other endpoint of x^ as t^. From t^ choose 
y^ different from existing links in G such that g^>0 and designate the 
other endpoint of y^ as t^. If no such y^ exists, go to Step 5(<i)» 

3. Let i = i + 1. Choose x^ (which currently joins t ^. ^ to 
t^^) and y^ as follows: 



42 



a. is chosen so that it is the longest link in G 

incident to vertex t^. , 

2i-l 

b. is some link that can be drawn from vertex t^^ sub- 
ject to (c), (d), and (e) below. If no y^ exists, go to Step 4. 

c. To guarantee that x's and y's axe disjoint, x^ cannot 

be a link previously joined (i.e,, ay., j < i) , and similarly y. cannot 

J 1 

be a link previously broken, 

i 

d. P. = S g. > 0, (Gain cilterion) 

1 ^ J 

e. y. chosen must permit breaking of an x. , > 

1 1+1 

f. Before y^ is drawn, check if closing up by joining 

t^^ to t^ will give a gain value better than the best obtained previously. 
Let y^* be a link connecting t^^ to t^ and let g^* = | ^i | " I ^i* | * 

If P. . + g. * > P*, set P* = P, ^ + g.* and let k = i. (P* is nov; the 
best improvement in G. P* > 0 and must be monotonically nondecreasing. 
The index k defines the sets to be exchanged to achieve P*) . 

4. Terminate the search for x^ and y^ in Steps 1 to 3 when 

either 

a. no further links x^ and y^ satisfy 3(c) to (e). 

b. P^ •< P* 

Delete the x^'s and replace them by the y^'s, i = l,...,k, 

to obtain G' , 

5. Test the feasibility of G' (see Sec. 4.2.3)* If feasible, 
take G' as the starting network and repeat from Step 1, Otherwise gen- 
erate another initial solution using the Starting Routine. 

6. If P* = 0, a limited backtracking is invoked; 

a. Repeat Steps 3 and 4, choosing y 2 's in increasing cost 
as long as they satisfy the gain criterion g^ + g 2 > 0* 
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b. If all choices of in Step 3(b) are exhausted vdthout 
profit, return to Step 3(a) and try the alternate choice for x^o 

c. If this also fails to give improvement, a further back- 
up is performed to Step 2, where the y^^'s are examined in increasing cost. 

d. If the y^^'s are also exhausted without profit, try the 
alternate x^ in Step 1. 

e. If this fails, select another vertex t^ and repeat 

from Step 1. 

7. The routine terminates when all vertices in G have been 
designated t^ vdthout profit. Return to the Starting Routine. 

A justification of Step 6 is in order; Experimental results 
by Lin and Kernighan showed that for the 2-redundancy network, backtrack- 
ing on the third level and higher, i.e., i > 3» indicated considerable 
time penalty, hence the only alternates for x^ , y^ and ^£'^2 first 

and second levels are examined. Other refinements to the procedure dis- 
cussed in the original paper are not applied here since actual computer 
runs may indicate different time-cost trade-offs for networks with redun- 
dancy greater than 2. 

4,2,3 Feasibility Test 

The optimizing routine for 2-redundancy network has a built-in 
feasibility test but is not applicable to higher redundancy networks. 

The fact that G is feasible does not imply that G' is also feasible even 
though the degree of the vertices satisfies the req.uirement. It is there- 
fore necessary to test the redundancy after every transformation. Rather 
than applying Algorithm 2 to examine the new network G' , the number of 
flow calculations can be reduced significantly using the theorems stated 
by Steiglitz, etal for an X-change or 2-change. Extending this to a k- 
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change transfoimation and vdth uniform redundancy, we need only to check 
the redundancy between the pairs of vertices whose connecting links were 
deleted, AlgoritVim 1 is most suitable for investigating feasibility after 
every any transfoirmation. 

Example ; 

Given a set of stations whose geographical arrangement is trans- 
ferred on an Euclidian plane as in Fig, We attempt to construct a 

minimum cost flow network vdth a uniform redundancy of 3* We model the 
network as a graph of vertices and branches and illustrate the applica- 
tion of the design procedure. 

Start by numbering the vertices uniformly at random as in Fig, 
Measure the distances between vertices (assuming proportional, to cost) and 
obtain the cost matrix. Add branches one at a time and keep track of the 
number of incident branches on each vertex by using an updated require- 
ment table as shown. Row 1 of this table is filled up by the assigned 
vertex requirements. This row has all the same entries so start with 
vertex 1, The cost matrix shows branch (1,6) has the least cost (here 
the parentheses are used in place of brackets to represent undirected 
branches) , Encircle this entry in the matrix and construct the branch 
(1,6) in Fig, 5*2, Form row 2 of the requirement table by reducing en- 
tries for vertices 1 and 6 by unity. The highest updated requirement in 
this row is 3 "the highest on the list is vertex 2, From the cost ma- 

trix we see that branch (2,7) has the least cost. Encircle this entry 
and construct the branch (2,7) t The procedure is continued and the se- 
quence of branches constructed next are (3,^), (5,8) ,(l,7) »(2,8) ,(3»6) , 
(^,5), (1,2), (3, 7), (4,8), ajid (5,6), 
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Fig. 5.3 



Fig. 5.4 
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When the number of vertices is odd, the last row of the requirement 
table has one vertex left with an updated requirement of 1 while the rest 
are zero. Connect this vertex to one which produces the least branch cost, 
at the same time satisfying Step 5(c) of the starting routine. 

The initial solution generated by the starting routine is shown by 
solid lines in Fig, 5*2. It is found feasible using Algorithm 2. The 
cost of the network is obtained by summing up the values of the encircled 
entries in the cost matrix. For this example, the initial cost is 261. 

We apply now the Optimizing Routine to reduce the cost. Set P* = 

0 and i = 1, In Fig. 5*2, choose vertex 1 as t^. x^ is selected such 
that it is the longest among the three incident branches to t^. Vertex 

2 which is the other endpoint of x^ is designated t^. From t^ find a y^ 

using Step 3 of the routine such that g^^ > 0. Only the imaginary branch 
(2,4) is available since the other possibilities (2,5)j(2,6), and (2,3) 
would result to g^ < 0. Choosing (2,4) as y^ the gain = I 1 “ 

|y^j = 3« The overall profit = 3 > ?*• Designate vertex 4 as ty 

Now let 1=2, Choose x^ as indicated. Assign t^ to vertex 8 

which is the endpoint of x^. Using Step 3 again, is selected to be 
the imaginary branch (8,?). The gain §2 = 9 and ^2 ^2 ~ 

fore going further, check if joining t^^ to t^ would improve P*, With 
( 8 , 1 ) as y 2 *» §2* ~ ^ P* = Pj^ + g 2 * = 12. We have to terminate the 
, search for x^'s and y^'s since P* = P 2 (Step 4). Had P* been greater 
than P 2 we could have continued by designating vertex 7 as t^ and look 
for Xj and y^. 

Construct the Y's and delete the X' s. The network with improved 
cost of 249 (not yet optimal) is shown in Fig. 5»3« The cost is obtained 
by subtracting the costs of x^ and X 2 from the initial value and then 
adding the costs of the branches y^ and 
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Feasibility of the new network is tested using Algorithm 1 since 
we need only to check the branch connectivity between the vertex-pairs 
1-2 and ^-8. This network is feasible so we can iterate the procedure 
using Fig, 5*3 as the new starting solution. 

Satisfying the vertex requirement does not necessarily mean fea- 
sibility is preserved. Fig. 5*^ obtained by the optimizing procedure is 
not feasible as it has a minimum connectivity of 2. 
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V. CONCLUSION 



Analysis procedures provide convenient tools for evaluating the 
"strength" or "weakness" of flow networks with respect to some physically 
meaningful criteria. The computer analysis programs were prepared as im- 
plementation of the algorithms developed and are capable of investigating 
large networks based on either branch or vertex survivability measure. 

It was shown that the analysis tools are not only used for measuring 
performance of existing systems but are also applied effectively in their 
design particularly when either explicit solutions are impractical or so- 
lutions do not appear to be obtainable. A heuristic algorithm for the 
design of minimum-cost survivable networks is presented. Analysis tech- 
niques are found indispensable in the procedure. 

Further study is indicated in design perhaps by encoding the algorithm 
for computer running and to investigate the folloi-ang: 

1. V/hat level of backtracking vjould be considered optimal? 

2. How often does infeasibility occur in the starting routine? The 
probability of occurrence may be different from that obtained by Stelglitz 
since the criteria are not the same. 

3. ' Are refinements described in [l6l applicable to networks with 
redundancy greater than 2? 
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APPENDIX A 



FLOWCHARTS 



A-1. s-t brainIch connectivity 
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A-2. MULTITERMINAL BRANCH CONNECTIVITY 



( Undirected Network ) 




52 



A- 2. (continued) 
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A-3. MULTITERMINAL BRANCH CONNECTIVITY 
( Directed Network ) 




54 



A-3. (continued) 
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A-4, s-t VERTEX CONNECTIVITY 
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• no 



Search for 


vertex 
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A-5. UNIFORM R- CONNECTIVITY 



© 



no 










X 

Read, N, R, 
connection matrix 






) ^ 




Choose source 
vertex ■ 








yes 



Choose a sink vertex 

V.. List no. of ver- 
3 

tices previously des- 
ignated sink that are 
adjacent to v 









Network is 
R-connected 



Remove all the pre- 
viously designated 
sink vertices adja- 
cent to V . and verify 
0 

R-connectivity be- 
tween V and V . using 
s 3 

Algorithm 4. 




no 



57 



ooo oooooooooooooooooooooooooooooooooooooooooooonoo 



APPENDIX B 



COMPUTER PROGRAMS 



B-1. PROGRAM HAXCON-1 



THIS PROGRAM COMPUTES FOR THE MULTITERMINAL ERANCH- 
CCNNECTIVITY OF A DIRECTED OR MIXED NETWORK. THE NETWORK 
IS REPRESENTED BY ITS CONNECTION MATRIX FOR COMPUTATION. 
THE NETWORK MUST HAVE CNLY SINGLE BRANCH BETWEEN VERTICES. 

THE PROGRAM PROCEEDS BY READING-IN THE INPUT CATA CON- 
SISTING OF THE TOTAL NUMBER OF VERTICES AND THE CONNECT- 
ION MATRIX. IT THEN DETERMINES THE INITIAL S’-T CUT WITH 
THE FIRST AND NTH VERTICES AS SOURCE AND SINK RESP. 

HAVING FOUND AN S'-T CUT, IT THEN FINDS OTHER S-T CUTS 
KEEPING THE SINK VERTEX CONSTANT AND VARYING ONLY THE 
SOURCE VERTEX FROM 2 TC N. AFTER COMPLETING THIS LOOP, 

THE SINK VERTEX IS SET TO N- 1 AND THEN FINOS AN S'-T CUT 
BETWEEN VERTICES 1 AND N-L. THE SOURCE IS AGAIN VARIED 
FROM 2 TO N-1 FDR THE OTHER S-T CUTS, 

THE PROCESS IS REPEATED UNTIL N(N-l) CUTS ARE DETER- 
MINED. EACH REPETITION INVOLVES SUBROUTINE FLCVAR. 

THE OUTPUT IS A PRINTOUT OF THE NXN BRANCH-CONNECTIVITY 
MATRIX. 

IDENTIFIER NOMENCLATURE 



ARRAYS 

KMAT— INPUT CONNECTION MATRIX WITH ENTRIES 0, 1,-1 

OMAT S»-T FLOW MATRIX REPRESENTING AN S'-T FLOW 

PATTERN. THIS MATRIX IS UPDATED BEFORE A NEW SET 
OF S AND S' VERTICES ARE DESIGNATED. 

FLOMAT- S-T FLOW MATRIX REPRESENTING S-T FLOW PAT- 
TERN RESULTING FROM AUGMENTATION AND HOMING STEPS. 
NONZERO ENTRIES OF THIS MATRIX UPDATE CMAT. 

TERMAT- MATRIX DENOTING THE MAXIMUM BRANCH-CCNNEC T- 
IVITY BETWEEN ANY VERTEX-PAIR. 



VECTOR 

KCDE- 



— LIST OF VERTICES OF THE NETWORK GRAPH, 



ADDITIONAL INFORMATION 



THIS PROGRAM IS COOED IN FORTRAN LANGUAGE. NETWORK OF 
ANY SIZE CAN BE HANDLED BY THIS PROGRAM WITHIN THE LIMIT- 
ATIONS OF THE COMPUTER STORAGE. THE DIMENSION, READ, AND 
WRITE STATEMENTS MUST HOWEVER BE ADJUSTED ACCORDINGLY. 



MAIN PROGRAM 



DIMENSION KMAT(25,25), FLOMATI 25 ,251 , DMAT(25,25) 
DIMENSION TERMAT ( 25,25) 

DIMENSION NODE (25) 

COMMON /ARRl/ KMA T , FLOMAT , DMAT 
CCMMON /ARR2/ TE^(MAT 

COMMON /PT/ SOURCE , SPRIME , SI NK,AMARK ,BRCUT , N ,G 
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READ (5,9) N 

READ (5,10) ((KMAT(I,J) ,J=1 ,25) ,I = 1,N) 
INITIALIZE TERMINAL MATRIX 

DO 2 X=1,N 

DC i Y=l,25 
TERMAT( X,Y)=0 

1 CONTINUE 

2 CONTINUE 



COMPUTE FOR THE N(N~1) BRANCH CUTS 

DC 8 K=1,N 
J=N+1-K 

SET TO ZERO THE FLOW MATRICES 

DO 4 L=1,N 

DO 3 K=1,N 
CMAT(L,M )=0 
FLCMAT( L,M)=0 

3 CONTINUE 

4 CONTINUE 

SPRIME=0 

DO 7 1=1, N 

IF (I .EQ.J) GC TO 7 

SCURCE=I 

SINK=J 

BRCUT=0 

CALL FLCVAR 

ENTER CUT VALUE INTO THE TERMINAL MATRIX 
TERMATd , J)=BRCUT 
SPRIME=I 

TRANSFER FLOMAT ENTRIES INTO DMAT 
DO 6 L=1,N 
DO 5 M=1,N 

IF (FLCMAT(L,M) .EQ.O) GO TO 5 
DMAT(L,M)=FLOMAT(L,M) 

RESET FLOMAT TO ZERO 
FL0MAT( L,M)=0 

5 CONTINUE 

6 CONTINUE 

7 CONTINUE 

8 CONTINUE 



9 

10 

11 

12 



COMPUTATION 
WRITE (6,11) 
WRITE (6,12) 
CALL PRTMAT 
STOP 



COMPLETED 



FORMAT 

FORMAT 

FORMAT 

FORMAT 

END 



(16) 

(2513) 

( • O' ,T14, ‘ BRANCH-CONNECTIVITY 
(T14, ' 



MATRIX • ) 
• ) 
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c 

c 

c 

c 

c 

c 



c 

c 

c 



c 

c 

c 



c 

c 

c 

c 



c 



c 

c 

c 



c 



c 



c 

c 

c 



c 

c 



SUBROUTINE FLOVAR 

THIS SUBROUTINE DETERMINES AN S-T CUT ONCE AN INITIAL 
S'-T CUT IS FOUND. IT MAXIMIZES AN S~T FLOW USING THE 
EXISTING S»-T FLOW PATTERN. 

SUBROUTINE FLOVAR 
IMPLICIT INTEGER(A-Z) 

DIMENSION KMAT(25,25), FLOMAT ( 25 , 25 ) , DMAT(25,25 ) 
DIMENSION NUDE(25) 

COMMON /ARRl/ KMAT , FL OM AT » DMAT 
COMMON /CH/ NODE 

COMMON /PT/ SOURCE, SPRIMEtSINK, AMARK,BRCUT»NtG 
ERASE ALL VERTEX LABELS 

1 DC 2 J=1,N 
NODE! JI=0 

2 CONTINUE 



SEARCH FOR AUGMENTATION PATH 
CALL AUGMNT 
IF (AMARK) 3,24,15 

AN S-G AUGMENTATION PATH IS FOUND 
SEARCH FOR HOMING PATH 

3 Q=G 

4 DO 5 J=1,N 

IF (DMAT(Q,J}.EQ.l) GO TO 6 

5 CONTINUE 

GO TO L 

6 NOOE( J ) -Q 
Q = J 

IF (Q.EO.GJ GO TO 19 
IF (J. Eg. SINK) GO TO 7 
GO TO 4 

A HOMING PATH IS FOUND 
SEARCH FOR TRUNCATION PATH 

7 Z = G 

IF (G. EQ. SPRI ME) GO TO 12 



8 DO 9 1=1, N , ^ 

IF (DMATd ,Z).EQ.l) GO TO 10 

9 CONTINUE 



10 



IF (Z.EQ.G) GO TO 12 

GO TO 11 

NCDE(1)=Z 



IF (Z.EQ.G) GO TO 21 
IF (Z .EQ.SPR IME ) GO TO 
GO TO 8 



11 



A TRUNCATION PATH IS FOUND 
TRUNCATE SPRIME-G PATH 

11 I=NODE(Z) 

DMAT(Z, I) = 0 
Z=I 

IF (Z.EQ.G) GO TO 12 
GO TO 11 

TRACE AND ERASE G-T PATH 

12 Q=SINK 

13 J=NODE(G) 

FLOMAT (J,Q)=1 
OMAT( J, Q)=0 

Q = J 

IF (J.EQ.G) GO TO 14 
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GO TO 13 

AUGMENT S~G PATH 

14 IF (G.EQ. SOURCE) GO TO 23 
J =G 

GO TO 16 

AUGMENT S-T PATH 

15 J=SINK 

16 I=NOOE(J) 

IF (I .LT.O) GO TO 17 

INCREASE FORWARD FLOW BY UNITY 
FLOMAT( I, J ) = 1 
GO TO 18 

^DECREASE BACKFLOW BY UNITY. 

FLCMATI J,I )=0 

18 J=I 

IF (J.EQ. SOURCE) GO TO 23 
GO TO 16 

A LOOP OF FLOW IS FOUND BY HOMING ROUTINE 

19 Q = G 

20 J=NODE(C) 

DMAT( J»Q) =0 
Q = J 

IF (O.EQ.G) GO TO 1 
GO TO 20 

A LOOP OF FLOW IS FOUND BY TRUNCATION ROUTINE 

21 Z=G 

22 I=NODE(Z) 

DMATIZ, I )=0 
Z = I 

IF (Z.EQ.G) GO TO 1 
GC TO 22 

23 BRCUT=BRCUT+1 
GO TO 1 

NO MORE AUGMENTATION PATH IS FOUND. FLCW PATTERN 
IS MAX 

24 RETURN 
END 



SUBROUTINE AUGMNT 

THIS SUBROUTINE IS USED TO DETECT AN S-G OR S-T PATH 
USING THE MODIFIED FORD AND FULKERSON LABELING ALGORITHM 

SUBROUTINE AUGMNT 
IMPLICIT INTEGER(A-Z) 

DIMENSION KMAT(25,25), FLOMAT ( 25 t 25 ) t 0MAT(25»25) 
DIMENSION NODE(251f JST0RE(10)i ISTORE(IO) 

CCMMON /ARRl/ KMAT , FLOMAT , DMAT 
COMMON /CH/ NODE 

CCMMON /PT/ SOURCE,SPRIMEtSINKT AMARK,BRCUT,N,G 

I=SOURCE 
L = 0 
M=0 
Z=0 

AMARK=0 

NCDE( I)=S0URCE . 

SEARCH FOR AUGMENTING PATH 
1 DO 4 J=1,N 
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c 



TC BRANCH BETWEEN VERTEX PAIRS 

I F (DMAT( J ,I ) .EQ.OI GO TO 2 
IF (J .EQ. SOURCE) GO TO 4 

Tti BACKFLOW IS ENCOUNTERED 

IF (NODE! J) .NE.O) GO TO 4 

LABEL VERTEX J 
NCD£( J) =-I 
GC TO 3 

IF (KMAT(I ,J) .NE. 1) GO TO 4 

A FORWARD BRANCH IS ENCOUNTERED. TEST FCR FORWARD 
FLOW 

IF (DMAT( I t J ) .NE.O ) GO TO 7 
IF (FLOMATI I , J) „NE.O) GO TO 4 
IF (NODE( J ) ,NE.O ) GO TO 4 

LABEL VERTEX J 
NODE(J)=I 



IS THE SINK VERTEX REACHED? 

3 IF (J.EQ.SINK) GO TO 8 
M = M+1 

JSTORE(M) =J 

4 CONTINUE 

IF (Z.NE.O) GO TO 6 
IF (M.EQ.O) GO TO 9 



TRANSFER CONTENTS OF JSTORE TO ISTORE 



DC 5 S=1,M 
ISTOREI S)=JSTORE( S) 

5 CONTINUE 

Z=N 
M = 0 

6 I=ISTORE(Z) 

1 = 1-1 

GO TO 1 

AN S-6 AUGMENTATION PATH IS FOUND 

7 G=I 
AMARK=-1 
GO TO 10 

AN S-T AUGMENTATION PATH IS FOUND 

8 A MARK =1 
GO TO 10 

NO AUGMENTATION PATH FOUND 

9 AMARK=0 
10 RETURN 

END 



SUBROUTINE PRTMAT 

THIS SUBROUTINE PRINTS TERMAT IN PROPER FORM AND IN 
FORMAT 3513. 

IF ARRAY LENGTH IS .LE. 35 A SINGLE LOOP SECTION PRINTS 
THE MATRIX, OTHERWISE IT PRINTS THE FIRST 35 COLUMNS FOL- 
LOWED BY THE NEXT 35 COLUMNS AND SO ON IN MULTIPLES OF 35 
USING MULTIPLE LOOPS. THE REMAINDER COLUMNS ARE PRINTED 
IMMEDIATELY BELOW THE LAST 35-COLUMN GROUP. 
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SUBROUTINE PRTMAT 
IMPLICIT INTEGER(A-Z) 

DIMENSION TERMAT (25 ,251 
COMMON /ARR2/ TERMAT 

^ CCMMCN /PT/ SOURCE, SPRIME, SINK, AHARK,BRCUT 

C HOW MANY MULTIPLES OF 35 IS N? 

MULPLE=N/35 

IF (MULPLE.GE .1) GO TO 2 

THIS SECTION PRINTS ARRAY WITH LENGTH 
DO 1 I=1,N 

V^RITE (6,81 I, (TERMAT(I,J), J=1,N) 

1 CONTINUE 

GO TO 7 

THIS SECTION PRINTS ARRAY WITH LENGTH 
IN GROUPS OF 35 

2 M = 1 
K=1 

3 MFLS35 = M*-35 
OO A I=1,N 

WRITE (6,8) I, (TERMAT ( i ,J 1 , J=M, MPLS35) 

A CONTINUE 

M = M-f-3 5 

IF (K-EQ.MULPLE ) GO TO 5 
WRITE (6,9) K 
WRITE (6,10) 

K = K-!-l 
GO TO 3 

THIS SECTION PRINTS REMAINDER COLUMNS 

5 REMDR=N-MULPLE^35 
IF (REMDR.LE.Ol GO TO 7 
WRITE (6,11) 

WRITE (6,12) 

C 

DO 6 L=1,N 

WRITE <6, 8) L, (TERMAT(L,J) , J=M,N) 

6 CONTINUE 
C 

7 RETURN 
C 

8 FORMAT (>0‘V T8, 12, •/', IX, 3513) 

9 FORMAT (//TIA, 'CONTINUATION ', 13, 2X, 'OF 

1C FORMAT ( TIA, ' ' ) 

11 FORMAT (//TIA, 'LAST SECTION UF TERMAT') 

12 FORMAT (TIA,' — •) 

END 



, N , G 



.LT. 35 



.GE. 35 



TERMAT' ) 
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PROGRAM MAXCOK'-IA 



C THIS PROGRAM COMPUTES FCR THE BRANCH- CONNECT I V ITY BET- 
C WEEN ANV TWO GIVEN VERTICES OF A DIRECTED OR MIXED NET- 
C WCRK. THE NETWORK GRAPH IS REPRESENTED BY ITS CONNECTION 
C MATRIX FOR CGRiPUTATI ON. NETWORK MUST HAVE ONLY SINGLE 
C BRANCH BETWEEN VERTICES. 

C 

C THE NUMBER OF VERTICES, THE SELECTEC SOURCE AND SINK 
C VERTICES, AND THE CONNECTION MATRIX ARE THE INPUTS. THE 
C PROGRAM; DOES A MAXELOW CALCULATION USING THE MOCIFIED 
C FORD AND FULKERSON LABELING ALGORITHM. COMPUTED CUT VALUE 
C IS THE NETWORK CONNECTIVITY WITH RESPECT TO THE SPECIFIED 
C VERTEX-PAIR, 

C 

C SEVERAL VERTEX-PAIRS MAYBE LISTED WITH THE INPUT AND 
C THE PROGRAM DOES SIMILAR MAXFLOW CALCULATION FOR EACH 
C PAIR. 

C 

C THE OUTPUT IS A PRINTOUT CONSISTING OF THE GIVEN SOURCE 
C AND SINK, THE CONNECTIVITY VALUE, AND THE BRANCH-CUTSET 
C AS DENOTED BY THE VERTICES TO WHICH EACH BRANCH CUT IS 
C CONNECTED. 

C 

C THE PRESENT DIMENSION AND READ STATEMENTS MAYBE ADJUST- 
C ED ACCORDING TO THE SIZE OF THE NETWORK TO BE ANALYZED. 

C 

C THIS PROGRAM MAYBE PREFERRED OVER MAXCON-1 WHEN ONLY A 
C FEW SOURCE-SINK RAIRS ARE TO BE ANALYZED. 

C 

C 

IMPLICIT INTEGER(A-Z) 

DIMENSION KMAI( 50,50), FLOMAT 1 5 0 , 5 0 ) , NOCEISO), NX(50) 
DIMENSION N0TNX(20), JSTORE(20), ISTOREI 20 
C 

READ (5,221 N 

READ (5, 23] ( ( KM AT ( I , J ) , J = 1 , 50 ) , I = 1 , N ] 

WRITE (6,25) 

1 READ (5,2A) SOURCE, SINK 
IF (SOURCE. EQ.O) GO TO 21 

INITIALIZE FLOW MATRIX TO ZERO 

DO 3 I=1,N 

DC 2 J=1,N ■ 

FLOMAT( I , J) =0 

2 CONTINUE 

3 CONTINUE 



ERASE VERTEX LABELS 

A DC 5 J=1,N 
NODE( J] =0 
5 CONTINUE 



I=SOURCE 
L=0 
M = 0 
Z = C 

BRCUT=0 

NCDE( n =SOURCE • 

SEARCH FOR AUGMENTING PATH 

6 DO 10 J=1,N 
C 
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C TEST FOR CONNECTING BRANCH BETWEEN VERTEX-PAIR 

I F (KMAT( I,J ) ) 7, 10,8 

A BACKWARD BRANCH IS ENCOUNTERED. TEST FCR POSI- 
TIVE BACKFLOW. 

7 IF (FLOKATI J, I) . EC.O) GO TO 10 
IF (N0DE( J J .NE.O) GO TO 10 

LABEL NODE(J) 

N0DE( Ji =-I 
GO TO 9 

A FORWARD BRANCH IS ENCOUNTERED. TEST FCR ZERO 
FORWARD FLOW. 

8 IF (FLOKATd, J) .NE.O) GO TO 10 

IS N0DE<J) LABELED? 

IF (NODE( J) .NE.O ) GO TO 10 
NC0E( J) =I 

IS THIS THE SINK VERTEX? 

9 IF (J.EQ.SINK) GO TO 13 
H = M + 1 

J STORE ( iX) =J 

10 CONTINUE 

IF (Z.NF.O) GO TO 12 
IF (M.EQ.O) GO TO 16 

TRANSFER CONTENTS OF JSTCRE TO ISTQRE 

CO 11 S=1,M 
ISTORE( SJ=JSTCRE(S) 

11 CONTINUE 

Z = M 
f-dO 

12 I = ISTOPE(.Z) 

Z = Z-1 
GO TO 6 

WE REACH HERE IF AN AUGMENTING PATH IS FOUND 
ADJUST FLOW ON ALL BRANCHES ALONG TFIS PATH. 

13 I=NODE{J) 

IF < I .LT.O) GO TO 14 

INCREASE FORWARD FLOW BY UNITY 

flcmat( I , ji =1 

GO TO 15 

DECREASE BACKFLOW BY UNITY 

14 I=-I 
FLOMAT ( J, I ) = 0 

15 J = I 

IF (J .EO. SOURCE) GO TO 4 
GC TO 13 

THIS STEP IS REACHED IF LABELING AND FLCW AUGMENT- 
ATION IS COMPLETED. SEARCH FOR BRANCH CUTS. 

16 00 19 1=1, N 
IF (NODEd I.EQ.OI GO TO 19 

DO 18 J=1,N 
ACD = 1 

IF (KMAT( I , JKEQ.O) GC TO 18 
IF (NODE( J) .NE.OI GO TO 18 

IF FLCW OCCURS IN BOTH DIRECTIONS, NET FLCW IS ZERO 
IF (FLOHAT(I,JI .EQ.FLOMATI J,I) ) GO TO 18 
IF {FLOMATd, J) .EQ.IJ GO TO 17 
ADO=FLOMAT( I , J) 
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c 



17 



18 



BRCUT=BRCUT+ADO 



L=L+1 



N X { L) = I 
NCTNX (L)=J 
CGiNTINUE 



19 CCNTINUE 



COMPUTATION TERMINATES SINCE NO MORE AUGMENTING 
PATH IS FOUND. 

IF (BRCUT.NE.O) GO TO 20 
tvRITE (6t26J SOURCE , SINK, BRCUT 
GO TO 1 

20 WRITE (6,27) SOURC E , S I NK , BR CUT , ( N X { J ) , NOTN X ( J ) , J = 1 , L ) 
GO TO 1 

21 STOP 

22 FORMAT (16) 

23 FORMAT (2513, /25I3) 

2A FORMAT (216) 

25 FORMAT ( » C' , lOX, • SOURCE^ , lOX, ‘SINK*, lOX, • BR- 

1 CONNECTIVITY', lOX, 'BR-CUTSET*) 

26 FORMAT (‘O', 12X, 12, i3X , 12, 17X, 12) 

27 FORMAT {'.O', 12X, 12, 13X, 12, 17X, 12, 

1 (T67, 13, IX, • — 13/)) 

END 
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PRCGRAM MAXCON-2 



B-3. 



C THIS PROGRAM COMPUTES FOR THE MU LT ITEP. M I NAL 8RANCH- 
C CONNECTIVITY OF ANY GIVEN UNDIRECTED NETWORK. THE NETWORK 
C GRAPH IS PRESENTED FOR CUMPUTATICiN BY ITS CONNECTION MAT- 
C RIX. NETWORK CAN HAVE PARALLEL BRANCHES BETWEEN VERTICES. 

C 

C THE PROGRAM PROCEEDS BY FIRST PEADING-IN THE NUMBER OF 
C VERTICES, THE INITIAL VERTEX-PAIR CHOSEN ARBITRARILY, 

C and the CONNECTION MATRIX. THE CUT BETWEEN THE TWC 
C SELECTED VERTICES IS OBTAINED USING SUBROUTINE MINCUT, 

C THE CONNECTIVITY FOR THE CUT(X,XBARI IS THE CCMFUTED CUT 
C VALUE. THIS IS ENTERED INTO THE TERMINAL CAPACITY MATRIX 
C BY SUBROUTINE MATRXT. THE SUCCEEDING VERTEX-PAIRS ARE 
C SELECTED BY THE PROGRAM FRCM THE GENERALIZED X-VERTEX- 
C GRCUP. SUBRGUTiNc MINCUT IS EMPLOYED EACH TIME FOR 
C OBTAINING THE CUT VALUE. WHEN ALL VERTICES IN X ARE EX- 
C HAUSTED, THE X-VE RTF. X-GRDUP IS CONDENSED INTO A SINGLE 
C VERTEX BY SUBROUTINE MATRXC AND THE XBAF-VER TEX-GROUP IS 
C IN TURN PROCESSED SIMILARLY. FOR EACH. COMPUTATICN, SUB- 
C ROUTINE MATRXT UPDATES THE TERMINAL CAPACITY MATRIX. 

C 

C WITH N AS THE TOTAL NUMBER OF VERTICES, THE FPCGRAM 
C WILL CO N-1 SUCH COMPUTATIONS AND FINALLY PRESENT A MULTI- 
C terminal BRANCH-CCNNECTIV ITY MATRIX AS OUTPUT. 



C ....... . IDENTIFIER NOMENCLATURE 

C 

ARRAYS 

KMAT INPUT CONNECTION MATRIX WHICH HAS ENTRIES 

GREATER THAN UNITY INDICATING PARALLEL BRAN- 
CHES BcTWEEN VERTICES 

CHAT— j — MATRIX REPRESENTING THE NEW GRAPH FORMED RES- 
ULTING FROM VERTEX CONDENSATION. IT HAS ALSO 
ENTRIES GREATER THAN UNITY. 

FLOMAT — MATRIX INDICATING MAXIMUM FLOW THAT CAN EXIST 
BETWEEN ANY VERTlX-PAIR ALLOWING UNIT CAPA- 
CITY FDR EACH BRANCH. 

TERMAT — MATRIX DENCflNG THE MAXIMUM CONNECTIVITY BET- 
WEEN ANY VERTEX-PAIR. 

VRTX ) 

NCTVTXI -ARRAY STORAGE FOR SEQUENCE OF CUT S ( X , XBAR ) , 
HAVING MORE THAN ONE VERTEX ON EACH SIDE. 



C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



w p p "T p p r 

NX- --VECTOR STORAGE FOR X-VERTICES IN CUT(X,XBAR] 

CURRENTLY BEING PROCESSED. 

NGTNX VECTOR STORr^GE FOR X8AR-VERT ICES IN CUT 

(XrXBAR) CURRENTLY BEING PROCESSED. 

NLIST — -STORAGE FOR VERTICES WHICH HAVE OCCURED SING- 
ULAR IN A CUT{X,XBAR). 

CONSTANTS AND COUNTERS 

N TOTAL NUMBER OF VERTICES IN THE GRAPH 

SOURCE — CHCSEN SOURCE VERTEX 
SINK CHCSEN SINK VERTEX 

RCW -DESIGNATES ROW NUMBER IN ARRAYS VRTX AND 

NOTVTX 

BRCUT DENOTES CUT VALUE 

COUNT 1 — NUMBER OF VERTICES IN NX 
CCUNT2 — NUMBER OF VERTICES IN NOTNX 
NCOUNT — RUNNING TOTAL OF VERTICES IN NLIST 



C ADDITIONAL INFORMATION . 

^ THIS PROGRAM IS CODED IN FORTRAN LANGUAGE. THE PRESENT 
DIMENSION STATEMENTS REQUIRE I50K OF COMPUTER STORAGE. 
LARGER NETWORKS CAN BE HANDLED BY AOJliSTING THE DIMENSION, 
READ, AND WRITE STATEMENTS. 
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MAl.'J PROGRAM 

IMPLICIT INTEGLR(A-Z) 

DIMENSION KKAT( 35,35), FLQMAK 3 5 , 35) , CHAT(35,35) 
DIMENSION TFRMAT{35,35) 

DIMENSION VRTX(35 ,35I* NOTVTX (35 ,35) t N0DE(35) 
DIMENSION NX(35), N0TNX(35) 

DIMEInISION NLIST(35) 

COMMON /ARRIV KMA T , F LOMAT ,C MA T , VRTX , NOT VTX 
CCMMON /ARR2/ TERM AT 
COMMON /CHI/ KODE,NLIST 
CCMMON /CH2/ NX,NOTNX 

CCMMCN /PTl/ N, SOURCE, SINK, ROW, MARK, FLAG, BRCUT 
CCMMON /PT2/ COUNTl ,C0L'RT2 , NCCUNT 

READ (5,29) N, SOURCE , S INK 

READ (5,50) ( (KMAT( I , JI ,J = 1 ,35.) ,I=1,N) 

INITIALIZE VARIABLES 

DO 3 I=1,N 

DC 2 K=L,10 

DO 1 J=1,N 

CMAT( I , J)=KMAT ( I , J ) 

VRTX( I , J) =0 
N0TVTX( 1 , J)=0 
NLIST( J)^0 

1 CONTINUE 

2 CCKTINUE 

3 CONTINUE 



DO 5 1=1,35 

DO 4 J=1 ,35 
TERMATI I , J)=0 

4 CONTINUE 

5 CONTINUE 

BRCUT =0 

NCCUNT=0 

FLAG=0 

MARK=0 

R0W=0 

F IND INITIAL S-T CUT 

6 CALL MINCUT 

LIST DOWN LABELED VERTICES UNDER NX AND UNLABELED 
ONES UNDER NCTNX 
CCUNT1=0 
CCUNT2=0 

DO 7 J=1,N 
NX(J)=0 
NCTNXt J)=0 

7 CONTINUE 



C 



DC 14 I=1,N 

TEST IF THIS VFRTEX IS ALREADY IN NLIST 
IF ( NCOUNT.EQ.O) GO TO 9 



CO 8 JJ=1,NC0UNT 
IF (NLI ST( JJ) .EQ.I ) 
8 CONTINUE 



GO TO 14 
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TEST IF IT BELONGS TO SOME CONDENSED NCTVTX 
9 ISTCRE=POW 

1C IF ( I STORE.EQ.O) GC TO 12 
DC 11 KK=1,N 

IF (NOTVTX(ISTORE,KK).EQ.n GO TC 14 

11 CONTINUE 

I STGRE=IST0kE-1 
GO TO 10 

12 IF (N0DE( I) . EQ.O) GO TO 13 
CCUNT1=C0UNT1+1 
NXICOUNTl ) = I 

GC TO 14 

13 CCUNT2 = CCUNT2-H 
NCTNX (C0UNT2 )=I 

14 CONTINUE 



UPDATE TERMINAL CAPACITY MATRIX ENTRIES 
CALL MATRXT 

IF (COUNTl.FQ.l) GC TC 18 
IF (C0UNT2.EQ»1 I GO TO 17 

BOTH NX AND NOT NX HAVE MORE THAN ONE VERTEX 
MARK=0 
RCW=ROW+l 

STORE NX AND NOTNX IN VRTX AND NOTVTX RESP 

DC 15 LL=1,C0UNT1 
VRTX(ROW,LU =NX( LL I 
15 CONTINUE 



DC 16 LL=1,C0UNT2 
NOTVTXI ROW ,LU =NCTi\X(LL) 
16 CONTINUE 



CONDENSE NOTVTX INTO SINGLE NODE 
CALL HATRXC 



SELECT TWO VERTICES IN NX AS SOURCE AND SINK 
S0URCE=NX{ 1) 

S INK=NX (COUNTl ) 

GC TO 6 

NX HAS MORE THAN ONE VERTEX AND NOTNX HAS ONE 

17 NCCUNT=NCCUNT+1 
NLISTI NCCUNT) =N0TNX(1J 
SCUPCE=NX ( 1) 

SINK=NX(CCUNT1I 

GC TO 6 

NX HAS SINGLE VERTEX AND NOTNX HAS MERE THAN ONE 

18 NC0UNT=NC0UNT+1 
NLISl ( NCCUNT) = NX ( 1 ) 

IF (C0UNT2.EQ.1 J GO TO 19 
SCURCE=NOTNX( 1 ) 

SINK=NOTNX (C0LNT2 ) 

GO TO 6 

BOTH NX AND NOTNX HAVE SINGLE VERTEX 



C 



19 DC 20 LL=1, NCCUNT 

IF (NLISTILL ) .£Q. NOTNX 
2C CCNTINUE 



1 ) ) 



GO TO 21 
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KCCUNT=NCCUNT+1 

NLI ST( NCCUNT)=N0TNX(1J 



21 



I F 



TEST IF ALL 
( NCOUNT.EQcN) 



VERTICES ARE IN NLIST 
GO TO 28 



THIS SECTION PROCESSES 
PREVIOUSLY CONDENSED 
CCUNT1=0 
CCUNT2=0 



THOSE NOTVTX GROUPS 



CO 23 MM=1,N 

IF (NOTVTXIRGWtMM) ,EQ.O ) GO TO 22 
CGUNT2=CCUNT2+1 

22 IF (VRTX ( ROW T MM I . EQ.O ) GO TO 23 
CCUNT1=CCUNT1+1 

23 CONTINUE 



CONDENSE VRTX INTO SINGLE VERTEX 
RELABEL VERiICES ACCORDING TO PRESENT CUT 



DO 24 J=1,N 
NCDE( J) =--0 
24 CCNTINUE 



DC 26 I=lfC0UNT2 
DO 25 J=1,N 

IF (J.NE.N0TVTX(R0N,I I ) GO IC 25 
M0Dc{ J 1=1 
GC TO 2o 

25 CCNTINUE 

26 CONTINUE 

S0URCE = N0TVTX<R0W, II 
SINK=N0TVTX{ RUN, CUUNT2) 

ERASE PRESENT VRTX AND NOTVTX ROWS 

00 27 KK=1,N 
VRTXIROWrKKI =0 
NCTVTX ( ROW ,KK)=0 

27 CONTINUE 

MARK=1 
FLAG=1 
RCW=R0W-'l 
CALL MATRXC 
GO TO 6 

COMPUTATION TERMINATES 

28 WRITE ( 6,311 
WRITE (6,32) 

CALL PRTMAT 
STOP 

29 FORMAT (316) 

3C format (3512) 

31 FORMAT ('O', T14, 'BRANCH-CONNECTIVITY MATRIX') 

32 FORMAT (T14,' * ) 

END 



SUBROUTINE MINCUT 

THIS SUBROUTINE COMPUTES FOR TflE MAXIMUM CONNECTIVITY 
BETWEEN TWO SELECTFD VERTEX-PAIR USING THE MCCIFIED FCRC 
ANC FULKERSON LABELING ALGORITHM. 
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SUtkOUTINE H INCUT 
IMPLICIT INTEGER(A-Z) 

DIMENSION KMAT(35i35) ♦ F LOMAT { 3 5 , 35 J , CMAT(35,35) 
DIMENSION VRTX(35,35), N3TVTX ( 35 ♦ 25) , N00E(35) 
DIMENSICN J ST ORE (20), I ST ORE (20) 

DIMENSION NLIST(35) 

COMMON /ARRl/ KM AT , FLOM AT , CMAT , VR TX , NOT VT X 
COMMON /CHI/ NODEtNLIST 

COMMON /PTl/ N, SOURCE, SINK, ROW, MARK, FLAG, BRCUT 
COMMON /PT2/ C0UNT1,C0UNT2,NC0UNT 

INITIALIZE FLOW MATRIX TO ZERO 



DO 2 1=1, N 

00 I J=1,N 
FLOMATC I , J)=0 

1 CONTINUE 

2 CONTINUE 



ERASE VERTEX LABELS 



3 DC 4 J=1,N 
NODE! J) =0 

4 CONTINUE 



I = SOURC E 
M=0 
1 = 0 

BRCUT=0 

NCCE( I )=S0URCE 

SEARCH FOR AUGMENTING PATH 

5 DO 6 J=1,N 

TEST FOR CONNECTING BRANCH BETWEEN VERTEX-PAIR 
IF (CMAT( I ,J) .EO.O) GO TO 6 
IF ( FLOMATI I , J ) .EQ.CM AT( I , J ) ) GO TO 6 

IS NODE(J) LABELED? 

IF CNODEIJ ) .NE.O) GO TO 6 
NCDE( J) =I 

IS THIS THE SINK VERTEX? 

I F ( J.EQ. SINK) GO TO 9 
M=M-«-l 

JSTCRE(M)=J 

6 CONTINUE 

IF (Z.NE.O) GO TO 8 
IF (M«EQ.O) GO TO 10 

TRANSFER CONTENTS OF JSTORE TO ISTORE 

DO 7 S=1,M 
ISTQRE(S)=JSTCRE(S) 

7 CONTINUE 



Z = M 
M=0 

8 I=IST0RE(Z) 
Z=Z-1 
GO TO 5 



W£ PEACH HERE IF AN AUGMENlING PATH IS FGUNO 
ADJUST FLOW ON ALL BRANCHES ALONG 7FIS PATH 

^ FLOMAtI i'! J)=FLOMAT( I , J)+l 
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J = I 

IF (J,EQ. SOUPC?) GO TO 3 
GO TO 9 



THIS STFP IS RFACH5D IF LABELING AND FLOW AUGMENT- 
ATION IS S'^MPLETED. SEARfH FOR 6PANCH CUTS. 

10 DO 13 I =1 , N 

I F (NODE! I ) .FO.O) GO TO 13 

DO 12 J=1,N 
ADD = 1 

IF (CMAT( I , J ) .0) GO TO 1 2 

IF (NODE! J ) . NF.O ) GO th 12 

IF FLOW OCrjRS IN BOTH DIt'ECTIONS, MFT FLOW IS ZERO 
IF (FLOMA T( 1 , J) . cq.flomaT ( J , I ) ) GO TO 12 
IF (FL'^Mf.TU , J) .FQ. 1) GO TO 11- 
ADD =FLO MAT ( 1 , J) 

11 BRF|Jt=brcmit+ADD 

12 CONTINUE 

13 CONTI NUF 



RETURN 

FND 



. . . . SUBROUTINE MATRXT 

THIS SUBROUTINE ENTERS CO''^dJT':'D CU'^ VALU‘D INTO A Tc^RVI- 
NAL •'APACITY MATRIX tprmat /> mo FURTHER UPO-'T^c qthER EN- 
TRIES R'^SEP HM V^vtm'jm FI nw ALONG A GIVEN PATH. MATRIX 
T=FMAT BECOMES THE OjTPUf MUL t I t erm ] m AL CONNECTIVITY MAT- 
RIX ARTFR N-1 COMPUTATI hns. 

SUBPOUT If'R ^'ATRXT 
IMPLICIT IMTcgfr(A-Z) 

DIMPMSmu KMAT( 35,35), FL OMA T ( 3 5 , 35) , CMAT(35,35) 

D IMP NS! ON TRPMAT (35 ,3 5 ) 

DIMENSION VRTX( 35,tb) , NOT V TX ( 3 5 , 3 5 ) , N0DF(35) 
DIMENSION NX (35), M0TNX(35) 

DIMFNSION NLIST(35) 

COMMON /ARR]/ KMAT,FLOMAT,CMAT, VPTX,NOTVTX 
COMMP'N /ARR?/ l^RMAT 
COMMON /CHI/ NOOE,NLIST 
COMMON /RH2/ NXtNOTMX 

COMMON /PTl/ N, SNURC-P, S INK , ROW , MARK, FLAG, BRCUT 
COMMON /PT2/ COUNTl ,:0UNT2 , NfOUNT 

ENTER CUT VALUE 

DO 17 1=1, COUNTl 

00 16 J = 1,C.P!UNT2 
KUT1=NX ( I ) 

KUT2=N0TNX( J) 

TFPMAT ( KUT 1, KUT? ) =BRCUt 
tERMAT ( KUT2, K’JTl ) = BPCUT 

COMPARE CUT VALUE WITH OTHER ENTRIES 

00 15 LL=1,N 
IF (COUNTl .PQ.l ) GO TO 7 
IF (C0LINT2.F0.1 ) GO TQ 9 
IF (NODE( LL) .NF. 0) GO TO 15 
IF (FLAG.FQ.O) GO TO 8 

CONOEMSPD VRTX IS ON THE SOURCE SIDE OF THE CUT 

DO 1 KK = 1,NC0UN'P 

IF (iNLI ST( KK) .EQ. LL) GO TO 2 
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1 CONTINUE 

GC TO 3 

2 X=NLIST(1) 

IF (NODE(X) .EQ.O ) GO TO 8 
GC TO 15 

CONDENSED NOTVTX IS ON THE SOURCE SIDE OF THE CUT 

3 IKctr— RU// 

^ IF (IKEEP.EQ.O) GO TO 8 
DO 5 KK=l5N 

IF (N0TVTX(1KEEP,KK).EQ.LL> GO TO 6 

5 CONTINUE 

IKEEP=IKEEP-1 
GC TO 4 

6 XBAR=NOTVTX( IKEEPtl } 

IF (NOUEIXBAR I.EQ.O) GO TO 8 
GC TO 15 

7 IF (NODE( LU .NE. 0) GO TO 14 

8 IF (TERNAT (KUTlvLL j .LE.TERMAT(KUT1,KUT2)) GO TO 15 
TERMATI KUTl ,LU =BRCUT 

TERMAK LL,KUT1J=BRCUT 
GC TO 15 

9 IF (NODE(LL) *Nt.O) GO TQ 13 

DISCARD THIS VERTEX IF COMPONENT OF SCME CONDENSED 
VERTICES 

DC 10 JJ=l,NCOUNT 

IF (NLI STC JJj .EQ. LU GC TO 13 

10 CONTINUE 

I KEEP=R0W 

11 IF (IKEEP.EQ.O) GO TO 15 
DC 12 JJ=1,N 

IF (NOTVTX( IKEEP , JJ) .EQ.LL) GO TO 13 

12 CONTINUE 

IKEEP=1 KEEP-1 

IF ( IKEEP.EQ.O) GO TO 15 

GO TO 11 

13 IF (TERMAT(KUT2»LL).LE.TERMAT(KUT1,KUT2) ) GC TO 15 
TERMAT ( KUT2tLU = BRCUT 

TERMAT( LL rKUT2) =BRCUT 
GC TO 15 

14 IF {C0UNT2.EQ.1) GO TO 13 

15 CONTINUE 

16 CONTINUE 

17 CONTINUE 

RETURN 

END 



. . . SUBROUTINE MATRXC 

THIS SUBROUTINE CONDENSES A GROUP OF VERTICES ON ONE 
SICE OF A CUT INTO A GENERALIZED SINGLE VERTEX AND OB- 
TAINS A CAPACITY MATRIX CMAT REPRESENTING THE NEt^ GRAPH 
FORMED. 



SUBROUTINE MATRXC 
IMPLICIT INTEGER'! A-2) 

DIMENSION KMAT(35,35), FL OH AT ( 3 5 , 3 5) , CHAT! 35,35) 
DIMENSION VRTX(35,35), NOT VTX ( 3 5 , 3 5 ) , NODE (35) 
DIMENSION NLIST(35) 

COMMON /ARRl/ KMAT , FLOMAT , CMAT, VRTX, NOT VT X 
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COMMON 

COMMON 

COMMON 



/CHI/ NODEtNLIST 

-(PI N , SOURCE , S I NK , ROW , MARK » F LAG , BRC UT 
/PT2/ COUNT 1 ,CCUM2, KCOJKT 



INITIALIZE CAPACITY MATRIX TO ZERO 



DO 2 I = 1,N 

DC 1 J=1,N 
CMAT( I , Jl=0 

1 CONTINUE 

2 CONTINUE 



FILL UP CAPACITY MATRIX 
DO 13 I=1,N 

IF (NODt{ I J .EQ. 01 GO TO 13 
IF (NCGUNT.EQ.Oi GO TO 4 

CO 3 LL=l,NCOUNT 
IF (I ,EO.NLIST(LL) ) GO TO 13 

3 CONTINUE 

4 MELAG=0 
JFLAG=0 

DC 12 J=lvN 

IF (KMAT( I ,J) .EQ.OI GO TO 12 
IF (NODE! J J . EQ.O) GO TO 5 
CMAT( I , JJ,= KMAT( I , J j 
GO TO 12 

5 1NCEX=RCW 

I F (MARK.EQ. 1 ) GO TO 6 
IF (FLAG.NE.OJ GO TO 6 
GO TO 11 

6 DO 7 MM=1,NC0UNT 
IF ( J.EQ.NLIST(MM) ) GO TO 8 

7 CONTINUE 

GC TO 9 

PROCESSED VERTICES (NLISTI TO BE CONDENSED 

8 MFLAG = MFLAGH 
JTEMP=NLIS1 m 
CMAT( I , JTEMPJ=MFLAG 
CPAT( JTEMP,U =MFLAG 
GO TO 12 

UNLABELED VERTICES TO BE CONCENSED 

9 DO 10 KK=1,N 

IF (J.EQ.NOTVTXI INDEX ,KK) ) GC TO 11 

10 CONTINUE 
C 

INDEX=INDEX-1 

IF ( INDEX. EQ.O) GO TO 12 

GO TO 9 

11 JFLAG=JFLAG+1 
JTEMP=NOTVTX( INDEX, 1) 

CMAT( I , JTEMPI =JFlAG 
CMAT( JTEMPtl ) =JFLAG 

12 CONTINUE 
C 

13 CONTINUE 
C 

RETURN 

END 
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SUBROUTINE PkTMAT 

THIS SUBROUTINE PRINTS TERMAT IN PROPER FGRk AND IN 
FORMAT 3513. 



IF ARRAY LENGTH IS ,LE. 35 A SINGLE LOOP SECTION PRINTS 
THE MATRIX, OTHERWISE IT PRINTS THE FIRST 35 COLUMNS FOL- 
LOWED BY THE NEXT 35 COLUMNS AND SO ON IN MULTIPLES OF 35 
USING MULTIPLE LOOPS. THE REMAINDER COLUMNS ARE PRINTED 
IMMEDIATELY BELOW THE LAST 35-COLUMN GROUP. 



SUBROUTINE PRTMAT 
IMPLICIT INTEGERIA-Z) 

DIMENSION TERMAT (35,35) 

COMMON /ARR2/ TERMAT 

COMMON /PT),/ N, SOURCE, SINK, ROW, MARK, FLAG, BRCUT 

HOW MANY MULTIPLES OF 35 IS N? 

MULPLE=N/35 

IF (MULPLE.GE .1) GO TO 2 

THIS SECTION PRINTS ARRAY WITH LENGTH . LT . 35 



DC 1. I-1,N 

WRITE (6,8) I ,(TERMAT(I ,J) ,J=1 ,N) 

1 CONTINUE 

GO TO 7 

THIS SECTION PRINTS ARRAY WITH LENGTH .GE. 35 
IN GROUPS OF 35 

2 M=1 
K = 1 

3 MPLS35=M+35 



DC 4 I=1,N 

WRITE (6,8) I ,( TERMAT( I ,J) , J=M, MPLS35) 
A CONTINUE 



C 



C 

c 



M = M-}-35 

IF (K.EC.MULPLE) GO TO 5 
WRITE (6,9) K 
WRITE (6,10) 

K = K+1 
GO TO 3 

THIS SECTION PRINTS REMAINDER COLUMNS 

5 REMDR=N-MULPLE«35 

IF (REMDR.LE.O) GO TO 7 
WRITE (6,11) 

WRITE (6,12) 

DC 6 L=1,N 

WRITE (6,8) L,(TERMAT(LrJ) , J=M,N) 

6 CONTINUE 



7 RETURN 



8 FORMAT 

9 FORMAT 

10 FOP MAT 

11 FORMAT 

12 FORMAT 
END 



( ' 0* , T8 , 12 , • /' , IX, 

(//T14, 'CONTINUATION' , 



35 13) 
13, 2X, 



(T14, ' 

(/,/T14,'LAST SECTION OF TERMAT') 
(T14, • ' J 



OF TERMAT* ) 

• ) 
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PRCGF>AM MAXC0N-2A 



C THIS PROGRAM COMPUTES FOR THE MAXIMUM BRANCH- CCKNECT- 
C IVITY BETWEEN ANY TWO GIVEN VERTICcS OF AN UNDIRECTED 
C NETWORK* THE NETWORK GRAPH IS REPRESENTED BY ITS CCNNECT- 
C ION MATRIX FOR COMPUTATION. THE NETWORK CAN HAVE PARALLEL 
C BRANCHES BETWEEN VERTICES. 

C 

C THE PROGRAM STARTS BY READING-- IN THE NUMBER OF VERTICES, 
C THE SOURCE AND SINK VERTICES, AND THE CONNECTION MATRIX. 

C THE CONNECTIVITY OF A CUT(X,X3AR) IS DETERMINED USING THE 
C MODIFIED FORD AND FULKERSON LABELING ALGORITHM. 

C 

C THE OUTPUT CONSIST OF THE VALUE OF THE CONNECTIVITY AND 
C THE VERTEX PAIRS JOINING EACH BRANCH CUT. 

C 

C THIS PROGRAM MAYBE PREFERRED OVER MAXCON-2 WHEN ONLY 
C SEVERAL SOURCE-SINK PAIRS ARE TO BE INVESTIGATED. 

C 

C THE PROGRAM DIMENSION AND READ STATEMENTS MLST BE 
C ADJUSTED FOR LARGER NETWORKS. 

C 

C 

IMPLICIT INTEGER(A-Z) 

C DIMENSION KMAT(60,60), FLOMAT (60 , 60 ) , N0CE(60), NX(20) 

DIMENSION N0TNXI20I, JSTGRE(20J, ISTCRE(2CJ 
C 

READ (6,17) N 

READ (5,18) ((KMAT(I,J) ,J = 1,60) ,I=1,N) 

WRITE (6,20) 

1 READ (5,19) SOURCE, SINK 
IF (SOURCE. EQ.O) GO TO 16 

INITIALIZE FLOW MATRIX TO ZERO 

DC 3 I=1,N 

DO 2 J=1,N 
FLCMAT ( I, JJ-0 

2 CONTINUE 

3 CONTINUE 



ERASE VERTEX LABELS 

A DC 5 J=1,N 
NODE( J) =0 
5 CONTINUE 



I=SOURCE 
L=0 
M = 0 
Z = C 

BRCUT=0 

NCDE( I } =SOURCE 

SEARCH FOR AUGMENTING PATH 



6 DO 7 J=1,N 



IF 
I F 



TEST FCR CONNECTING BRANCH BETWEEN VERTEX-PAIR 
(KMATd ,J) .EQ. 0) GO TO 7 ^ ^ 

(FLOMAT ( I , J ) . EQ.KM AT ( I , J ) ) GO TO 7 



IS NOD£(J) LABELED? 

IF (NODE(J ) .NE.O) GO TO 7 
NGDE( J) =I 
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c 



c 



7 



IS THIS THE SINK VERTEX? 
IF (J.EQ.SINK) GO TO 10 



J STORE ( =J 
CGNT INU E 



IF {/.NE..0) GO TO 9 
IF (M.EQ.O) GO TO 11 



TRANSFER CONTENTS OF JSTORE TO ISTORE 

DC 8 S=1,M 
I STORE! S) =JST0RE( S) 

8 CONTINUE 



Z=M 
M = 0 

I=IST0RE(Z) 
Z=Z-1 
GO TO 6 



WE REACH HERE IF AN AUGMENTING PATH IS FOUND 
ADJUST FLOW ON ALL BRANCHES ALONG THIS PATH. 

10 I=N0DE(J) 

F L CMA T { I , J ) =F LCMA T ( I , J ) +1 
J= I 

IF (J.EQ. SOURCE) GO TO 4 
GO TO 10 

THIS STEP IS REACHED IF LABELING AND FLOW AUGMENT- 
ATION IS COMPLETED. SEARCH FCR BRANCH CUTS. 

11 DC 14 I=1,N 

IF (NODE! I ) , EQ.O) GO TO 14 

DC 13 J=1,N 
ADD=1 

IF CKMAT! I , J ) .EQ .0) GO TO 13 
IF {NODE( J) .NE.O) GO TO 13 

If FLOW OCCURS IN BOTH DIRECTIONS, NET FLOW IS ZERO 
IF {FLOMAT! 7 , J) .EQ.FLCMAT ( J , I ) ) GO TO 13 
IF {FLCHATd ,Jj .EQ.l) GO TO 12 
ACD--FLOMAT (I , J ) 

12 BRCUT = BRCUT-t-ADD 
L = L + 1 

NX(L) =I 
NQTNX(L)=J • 

13 CONTINUE 

14 CONTINUE 



C 



COMPUTATION TERMINATES SINCE NO MORE AUGMENTING 
PATH IS FOUND. 

IF (BRCUT.KE.O) GO TO 15 
WRITE (6,21) SOURCE, SINK, BRCUT 
G C TO 1 

15 WRITE (6,22) SOURC E , SI NK , BRCUT , ( NX ( J ) , NOT NX ( J ) , J = 1 , L ) 
GO TO 1 

16 STOP 



17 FORMAT (16) 

18 FORMAT (35 12, /25I2) 

19 FORMAT (216) ^ 

20 FCkMAT iOX, 'SOURCE', lOX, J^SINK 

12 , 
12 ^ 



21 

22 



1 



1 



FORMAT 

FORMAT 



CONNECTIVITY* , lOX, ' BR-CUTSET' ) 



( *0 
CO' 



END 



(T67, 13 



12X, 
12X, 
. IX, 



I 2, 

12 , 



13X, 12, 
13X, 12, 

, 13/)) 



17X, 

17X, 



ICA, 'BR- 
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PRCGRAM KAXCCIM-3 



C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



THIS PROGRAK CC.VPUTES FCR THE MULT ITERMIiMAL VERTEX- 
CCNNECTIVITY OF A DIRECTED OR MIXED NETWORK. THE NETWORK 
IS REFPESENTEC BY ITS CONNECTION MATRIX FCR COMPUTATION. 
THE NETWORK MUST HAVE ONLY SINGLE BRANCH BETWEEN VERTICES 

THE PROGRAM PROCEEDS BY READING- IN THE INPUT DATA CON- 
SISTING OF THE TOTAL NUMBER OF VE'^TICES AND THE CGNNECT- 
ICN MATRIX. IT THEN DETERMINES THE INITIAL S'-T CUT WITH 
THE FIRST AND NTH VERTICES AS SOURCE AND SINK RESP. 

HAVING FOUND AN S‘-T CUT, IT THEN FINDS OTHER S-T CUTS 
KEEPING THE SINK VERTEX CONSTANT AND VARYING CNLY THE 
SCURCE VERTEX FROM 2 TO N. AFTER COMPLETING TFIS LOOP, 

TEE SINK VERTEX IS SET TO N- 1 AND THEN FINDS AN S'-T CUT 
BETWEEN VERTICES 1 AND N-1 . THE SOURCE IS AGAIN VARIED 
FROM 2 TO N-1 FOR THE OTHER S-T CUTS. 

THE PROCESS IS REPEATED UNTIL N(N-II CUTS ARE DETER- 
MINED. EACH REPETITION INVOLVES SUBROUTINE FLCVAR, 

THE OUTPUT IS A PRINTOUT CF THE NXN VERTEX-CONNECTIVITY 
MATRJ X. 



IDENTIFIER NOMENCLATURE 



1 



ARRAYS 

KMAT INPUT CONNECTION MATRIX WITH ENTRIES 0, 1, 

DMAT S'-T FLOW MATRIX REPRESENTING AN S'-T FLOW 

PATTERN. THIS MATRIX IS UPDATED BEFORE A NEW SET 
OF S AND S* VERTICES ARE DESIGNATED. 

FLOMAT- S-T FLOW MATRIX REPRESENTING S-T FLOW PAT- 
TERN RESULTING FROM AUGMENTATION AND HOMING STEPS, 
NONZERO tNTRIES OF THIS MATRIX UPDATE CMAT . 

TERMAT- MATRIX DENOTING THE MAXIMUM BRANCH-CCNNECT- 
IVITY BETWEEN ANY VERTEX-PAIR, 

VECTOR 

NCDE LIST OF VERTICES OF THE NETWORK GRAPH. 



ADDITIONAL INFORMATION 



IHiS PROGRAM IS CODED IN FORTRAN LANGUAGE. NETWORK OF 
ANY SIZE CAN BE HANDLED BY THIS PROGRAM WITHIN THE LIMIT- 
ATICNS’OF THE COMPUTER STORAGE. THE DiMENSiOtN, READ, AND 
WRITE STATEMENTS MUST HOWEVER BE ADJUSTED ACCCRCINGLY, 



MAIN PROGRAM 



DIMENSION KiMAT(25,25)» F LOMAT ( 2 5 , 25 ) t DMAT(25,25I 
DIMENS ICN TERMAT ( 25,25) 

DIMENSICN NODE (25) 

COMMON /ARRl/ KMA T , F LOMA T , D MAT 
COMMON /ARR2/ TERMAT 

CCMMON /PT/ SCURCE , SPR IME , SINK, AMARK jVTXCLT ,N ,G ,F LAG 

READ (5 ,11) N , , K , 

READ (5,12) ( (KMATd, J) ,J = 1 ,I=1,N) 

INITIALIZE TERMINAL MATRIX 



DC 2 X=1,N 

DO 1 Y=l,25 
TERMAT (X,Y )=0 
1 CONTINUE 
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2 CONTINUE 



COMPUTE FOR THE N(N’-l) VERTEX CUTS 

DC 10 K=1,N 
J=N+1-K 



SET TO ZERO THE FLOW MATRICES 



DC 4 L=1,N 

DO 3 M=1,N 
CMAT(L,MJ=0 
FLCMAT { L, M) =0 

3 CONTINUE 

4 CONTINUE 

SPRIME^^O 

FLAG=0 

DC 9 I=lfN 

IF (I .EQ.J ) GO TO 9 

IS THERE A DIRECT CONNECTION BETWEEN SOURCE AND 
S INK? 

IF (KMAT( I , J) .GT.O I GO TO 5 

SOURCE^ 1 

S INK=J 

V1XCUT=0 

CALL FLOVAR 

ENTER CUT VALUE INTO THE TERMINAL MATRIX 
TERMATI I , JI=VTXCUT 
SPRIME= I 
60 TO 6 

5 TERMATI I , Jl =0 
GO TO 9 

TRANSFER FLOMAT ENTRIES INTO DMAT 

6 00 8 L=1,N 
DO 7 M=1,N 

IF (FLOMAT(L,M).EQ.O) GO TO 7 
DMAT(L, M)=FLCMAT (L,M) 

RESET FLOMAT TO ZERO 
FLOMAT ( L,M)=0 

7 CONTINUE 

8 CONTINUE 

FLA6=1 

9 CONTINUE 

10 CONTINUE 



COMPUTATION COMPLETED 
WRITE (6,13» 

WRITE (6,14) 

CALL PRTMAT 
STCP 
C 

11 FORMAT (16)^ 

1^ FCR(-9\T (^0^',T14 ,• VERTEX-CONNECT IV ITY MATRIX*) 

14 FORMAT IT 14,* 'J 

END 

C 
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.... SUBROUTINE FLOVAR . . 

THIS SUBROUTINE DETERMINES AN S-T CUT ONCE AN INITIAL 
S'-T CUT IS FOUND. IT MAXIMIZES AN S-T FLOW USING THE 
EXISTING S'-T FLCW PATTERN. 

SUBROUT INE FLOVAR 
IMPLICIT INTEGER! A-ZJ 

DIMENSION KMAT(25,25), FL OM AT ( 2 5 , 2 5 ) t DMAT(25,25) 
DIPENSICN NODE (25) 

CCMMGN /ARRl/ KMAT , FLOMAT , DMAT 
CCMMON /CH/ NODE 

CCMMON /PT/ SOURCE, SPR IME,S INK, AMARK,VTXCUT,N,G, FLAG 
ERASE ALL VERTEX LABELS 



1 DO 2 J=1,N 
NCDE( J) =0 

2 CONTINUE 



SEARCH FOR AUGMENTATION PATH 
CALL VTXCON 
IF (AMARK) 3,24,15 

AN S-G AUGMENTATION PATH IS FOUND 
SEARCH FOR HOMING PATH 

3 Q=G 

4 DC 5 J=1,N 

IF (DMAT(Q,J).EQ.l) GO TO 6 

5 CONTINUE 

GO TO 1 

6 NODE(J)=Q 
Q = J 

IF (Q.EQ.G) GO TO 19 
IF (J .EQ.SINK) GO TO' 7 
GC TO 4 

A HOMING PATH IS FOUND 
SEARCH FOR TRUNCATION PATH 

7 Z=G 

IF (G.EQ.SPR IME ) GO TO 12 

8 DC 9 I =1 tN 

IF (DMAT( I,Z) .EQ.l) GO TO 10 

9 CONTINUE 



1C 



IF (Z.EO.G) GO TO 12 
GC TO 11 
NODE! I ) =Z 



IF^ (Z.EQ.G) GC TO 21 
IF (Z.EQ. SPRIME) GO TO 
GO TO 8 



11 



A TRUNCATION PATH IS FOUND 
TRUNCATE SPRIME-G PATH 

11 I=NODE(Z) 

DMAT(Z, I)=0 

IF (Z.EQ.G) GO TO 12 
GO TO 11 

TRACE AND ERASE G-T PATH 

12 Q=S1NK 

13 J=NODE(Q) 

FLGMAT( J,Q)=1 
DMAT( J,Q)=0 

Q = J 

IF (J.EQ.G) GO TO 14 
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GO TO 13 

AUGMENT S'G PATH 

14 IE (G.EQ. SOURCE) GO TO 23 
J =G 

GO TO 16 

AUGMENT S-T PATH 

15 J=SINK 

16 i=NDDE(J) 

IF (1 ,LT.O) GO TO 17 

WEAKLY LABELEDe INCREMENT FRO FLOW BY 1 

F LCMAT ( I t J) =1 
GO TO 18 

17 STRONGLY LABELED.. DECREMENT BACKFLOW BY 1 

FLOMAT ( J, I ) = 0 

18 J = I 

IS SOURCE VERTEX REACHED? 

IF (J.FQ. SOURCE) GO TO 23 
GO TO 16 

A LOOP OF FLOW IS FOUND BY HOMING ROUTINE 

19 0 = G 

20 J=N0DEI0) 

DMAT( J, Q) =0 
Q = J 

IF (Q.EQ.G) GO TO 1 
GO TO 20 

A LOOP OF FLOW IS FOUND BY TRUNCATION ROUTINE 

2) Z=G 

22 I=N0DE(Z) 

DMATCZ, I)=0 
Z = I 

IF (Z.EQ.S) CO TO 1 
GC TO 22 

23 V1XCUT=V7XCUTH 
GO TO 1 

NO MORE AUGMENTATION PATH IS FOUND. FLCW PATTERN 
IS MAX 

24 RETURN 
END 



SUBROUTINE VTXCON ......... 

THIS SUBROUTINE LOCATES AN S-G OR S-T AUGMENTATION PATH 
USING FRISCH'S VERTEX-PAIR CONNECTIVITY ALGORITHM. 

IT IS INVOKED BY SUBROUTINE FLOVAR WHICH SPECIFIES THE 
SOURCE AND SINK VERTICES. THIS ROUTINE LABELS THE VERTI- 
CES WHILE LOCATING A PATH. INFORMATION WHETHER THERE IS 
AN AUGMENT! No PATH OR NOT IS RETURNED TO ROUTINE FLOVAR. 

SUBROUTINE VTXCON 
IMPLICIT INTEGER(A-Z) 

DIMENSION KMAT(25,25>t FLOMAT ( 25 , 2 5 ) t DMAT(25f25) 
OIMENSICN N0DE(25), JSTCRE(20)t 1ST0R£(20) 

COMMON /ARF-1/ KHA T , FL CM AT , DMAT 
CCMMON /CH/ NODE 

CCMMCN /PT/ SOURCE, SPR IME ,S INK, AMARK , VT XCUT ,N , G , F L AG 

INITIALIZE VARIABLES 
I = SOURCE 
L = 0 
M=0 
Z=0 
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AM4RK=0 

NCDE { I ) =SOURCE 



SEARCH FOR AUGMENTING PATH 



1 DO 6 J=1,N 

TEST FOR CONNECTING PATH RETl^JEEN VERTEX-PAIR 
IF (FLCMATI J,II .EQ.O) GO TO 3 
I F (NQOEI I i .LT.O J GO TO 2 

A BACKWARD BRANCH IS ENCOUNTERED. 

IF (J . EQ. SOURCE) GO TO 7 
IF (NODE! J ) .LT .0 ] GO TO 6 

VERTEXIJ) IS EITHER UNLABELED OR WEAKLY LABELED. 
STRONGLY LABEL THIS VERTEX. 

N0DE(J)==-I 
H=M + 1 

JSTORE ( M) =J 
GO TO 7 

2 IF (J .EQoSCURCE) GO TO 6 
NCDE( J) =-I 

GO TO 5 

3 IF (KMATU ,J) .NE.L ) GO TO 6 

A FORWARD FLOW IS ENCOUNTERED. TEST FOR ZERO FLOW. 
IF (FLAG.. Ey.O) GO TC 4 
IF (DMAT( I ,NE.O) GO TO 10 

4 IF (FLOMATUtJ) .NE.O) GO TO 6 
IF (NODL( JI.NE.O) GO TO 6 
NODE( J)=I 

IS SINK VERTEX REACHED? 

IF (J.EQ.SINK) GO TO 11 

5 M=M+l 

J STORE (M) 

6 CONTINUE 

7 IF (Z.NE.O) GO TO 9 
IF (M.EQ.O) GO TO 12 

TRANSFER CONTENTS OF JSTORE INTO ISTORE 

DO G S=lrM 
ISTOREI S) =JSTCRE (S) 

8 CONTINUE 

Z = f-' 

M = 0 

9 I=ISTORE(Z) 

Z = 2-l 

GO TO 1 

AM S-G AUGMENTATION PATH IS FCUNO 

10 G = I 
AMARK=-1 
GC TO 13 

AN S-T AUGMENTATION PATH IS FOUND 

11 AMARK=1 
GO TO 13 

NO AUGMENTATION PATH FOUND 

12 AMARK=0 

13 RETURN 
END 
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SUBROUTINE PRTMAT 



crouTf PRINTS TERMAT IN P-ROPER FCRP AND IN 

r L KMAT 35 1 3 • 



tpp^mat§tv ‘-OOP SECTION PRINTS 

PRINTS THE FIRST 35 CCLUMNS FOL- 
LCKED BY THE NEXT 35 CCLUMNS AND SO ON IN MULTIPLES OF 35 
USING MULTIPLE LOOPS. THE REMAINDER COLUMNS ARE PRINTED 
IMMEDIATELY BELOW THE LAST 35-COLUMN GROUP. 



SUBROUTINE PRTMAT 
IMPLICIT INTEGER(A-Z) 

DIMENSION TERMAT ( 25,25) 

CCMMON /ARR2/ TERMAT 

COMMON /PT/ SOUKC E, SPRIME , SI NK, AMARKtVTXCUT ,N ,G,FLAG 

HOW MANY multiples CF 35 IS N? 

MULPLE=N/35 

IF (MiULPLE.GE .1 ) GO TO 2 



THIS SECTION PRINTS ARRAY WITH LENGTH ,LT. 35 



DC 1 1=1, N 

WRITE (6,8) 1 , (TERMAT( I , J ) , J=1 ,N) 

1 CCNTINUE 

GO TO 7 

THIS SECTION PRINTS ARRAY WITH LENGTH .GE. 35 
IN GROUPS OF 35 

2 M = 1 
K = 1 

3 MlPLS35=M+35 
DC A I=1,N 

WRITE (6,3) I , (TERMAT( I ,J) , J = iM,MPLS35) 

4 CCNTINUE 

M = M-«-3 5 

IF (K.EQ. MULPLEl GC TO 5 
WRITE (6,9) K 
WRITE (6,10) 

K = K+1 
GC TO 3 



C 

C 

C 



5 



6 

7 

8 
9 

1C 

11 

12 



THIS SECTION PRINTS REMAINDER CCLUMNS 
REMDR=N-MULPLEa'-35 
IF (REMCR.LB.O) GO TO 7 
WRITE (6,11) 

WRITE (6,12) 

DO 6 L=1,N 

WRITE (6,8) L, (TERMAT(L, J) , J=M,N) 

CCNTINUE 

RETURN 

FORMAT ('O', T8, 12, '/', IX, 3513) 

FORMAT (//T14, ' CONTINUAT ION ' , 13, 2X, 'OF TERMAT*) 

FORMAT (T14,' — ') 

FORMAT (//T14,'LAST SECTION OF TERMAT') 

FORMAT (T14,' ') 

END 
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PRCGRAM MAXC0N-3A 



C THIS PROGRAM COMPUTES FOR THE VERTEX-CONNECTIVITY OF A 
C DIRECTED, AN UNDIRECTED, OR MIXED NETWORK USING FRISCH'S 
C VERTEX-PAIR CONNECTIVITY ALGORITHM. 

C 

C THE PROGRAM STARTS BY READING-IN THE INPUT CATA CONSIS- 
C TING OF THE NUMBER OF VERTICES, THE CONNECTION MATRIX, 

C AND THE SOURCE AND SINK VERTICES AS CHOSEN, THENCE IT 
C FINOS AN AUGMENTING PATH USING THE MODIFIED LABELING AL- 
C GCRITHM AND I NCRE AS ES / DECR E AS ES THE FLOW ONCE SUCH PATH 
C IS FOUND. THE PROCESS IS REPEATED UNTIL NO MG PE AUGMENT- 
C ING PATH IS FOUND. IT THEN SEARCHES FOR THE VERTEX CUTSET 
C BY FINDING A LABELED VERTEX AND AN UNLABELEO CNE WHICH 
C ARE BOTH CONNECTED AND HAVE UNITY F LCW BETWEEN THEM. 

C 

C THE OUTPUT IS A PRINTOUT OF THE SOURCE VERTEX, THE SINK 
C VERTEX, AND A LIST OF VERTICES COMPRISING THE CUTSET. 

C 

C FCR NETWOR^KS LARGER THAN THE DECLARED DIMENSION, THE 
C DIMENSION AND READ STATEMENTS MUST BE ADJUSTED. 

C 

C 

IMPLICIT INTEGER! A-Z) 

DIMENSION KMAT(50,50), FLOMATI 50 , 50) 

DIMENSION N00E(50J, JSTORE(20I, JSTORE{20J, VTXCUTUCJ 
READ { 5,241 N 

READ (5,25) ( (KMAT { I , J ) , J= 1 , 50) , I = 1 , N) 

WRITE (6,27) 

1 READ (5,26) SOURCE, SINK 
IF (SOURCE. EQ.O) GO TO 23 

IS THERE A DIRECT CONNECTION BETWEEN SOURCE 
AND SINK? 

IF (KMAT( SOURCE , SI NK) .GT.O) GO TO 22 

INITIALIZE FLOW MATRIX TU ZERO 

DO 3 I=-l,N 

DO 2 J=1,N 
FLOMAT ( I , J ) = 0 

2 CONTINUE 

3 CONTINUE 



ERASE VERTEX LABELS 

4 DO 5 J=1,N 
NODE( J) =0 

5 CONTINUE 



INITIALIZE VARIABLES 
I = SOURCE 
L = 0 
M=0 
Z=0 

NODEl I )=SOURCE 



SEARCH FOR AUGMENTING PATH 



6 DO 10 J=1,N 



IF 
1 F 



TEST FOR CONNECTING PATH BETWEEN VERTEX-PAIR 
(FLOf'lAT! J, I ) .EQ .0) GO TO 8 
(NODE! 1) .LT.O) GO TO 7 



A BACKWARD BRANCH IS ENCOUNTERED. 
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I F 
I F 



(J.EO. SOURCE) 
(NODL( J ) .LT.O) 



GO TO 11 
GO TO 10 



VEkTtX(J) IS EITHER UNLAbELED OR WEAKLY LABELED. 
STRONGLY LABEL THIS VERTEX. 

NODt(J) =-I 
M=M + 1 

JST0RE(M)=J 
GO TO II 

7 IF (J.EQ, SOURCE) GO TO 10 
NODE(J) =-I 
GO TO 9 



8 IF (KMATI I , J) .NE. 1 ) GO TO 10 



A FORWARD BRANCH IS ENCOUNTERED. 
IF (FLOMAT ( I , J ) .NE.O ) GO TO 10. 

IF (NODE! J ) .NE . 0) GO TO 10 
fxIODEI J ) = I 



TEST FCR 0 FLCW, 



IS SINK VERTEX REACHED? 

IF (J.EQ.SINK) GO TO 14 
9 M=M+1 

J STORE ( Mj =J 

10 CONTINUE 

11 IF (Z.NE.O) GO TO 13 
IF (M.EQ.O) GO TO 17 

TRANSFER CONTENTS OF JSTORE INTO ISTORE 

DC 12. S = UM 
I STORE I S) =JSTORE ( S) 

12 CONTINUE 



Z=M 
M = 0 

13 1=1ST0RE(Z) 
Z = Z-1 
GO TO 0 



WE REACH HERE IF AN AUGMENTING PATH IS FOUND. 

ADJUST ALL FLOWS .iLDNG THIS PATH. 

TEST LABEL ON VERTEX(J) 

14 I=N0DE(J) 

IF ( 1 .LT.O) GO TO 15 

VERTEX IS WEAKLY LABELED. INCREMENT FKC FLOW BY 1 
FLCMAT I It J) = l 
GO TO 16 

VERTEX IS STRONGLY LABELED. DECREMENT BACKFLOW 

15 I=-I 
FLOMATI Jr I )=0 

16 J=I 



IS SOURCE VERTEX REACHED? 
IF (J.EQ. SOURCE) GO TO 4 
GC TQ 14 



VERTEX LABELING AND FLOW AUGMENTATION IS COMPLTD. 
THIS SECTION SEARCHES FCR THE VERTEX CUTS 



17 DC 20 1=1 ,N 

IF (NODE! I ) .EQ.O) GO TO 20 
C 

DC 19 J=1,N 

IF (NOOE( J) .NE.O) GO TO 19 
IF (FLOMATd, J) .EO.O) GO TO 19 
IF (I .EQ. SOURCE) GO TO 18 
L=L+1 
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VTXCUK u =I 
GO TO 20 
IS L=L+1 

VTXCUK U =J 

19 CONTINUE 

20 CONTINUE 



COMPUTATION TERMINATES SINCE NO MORE AUGMENTING 
PATH IS FOUND. PRINT CUTSET. 

IF (L.EQ.O) GO TO 21 

ViRITE (6,28) SOURCE, SINK, (VTXCUTI J) ,J=1 ,L) 

GC TO 1 

21 WRITE (6,2 9) SOUPxCE,SINK 
GO TO 1 

22 WRITE (6,30) SOURCE, SINK 
GO TO 1 

23 STOP 

2''i FORMAT (16) 

25 FORMAT (2513/2513) 

26 FORMAT (216) 

27 FORMAT ('OS lOX, 'SOURCE', lOX, 'SINK', lOX, 

1 'VERTEX-CUTSET' ) 

28 FORMAT ('O', 12X, 12, 13X , 12, (T42, 413)/) 

29 FORMAT ( ‘ 0 ' , 12X, 1 2 , 13X, I 2 , T41, ' NO FEASIBLE FLOW') 

30 FORMAT ('O', 12X, 12, 13X, 12, T44, ‘O') 

END 
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PROGRAM MAXCON-^4 



C THIS PROGRAM VERIFIES THE R-C GNNE C T I VI T Y OF AN UNDIR- 
C ECTED NETWORK FOR A GIVEN UNIFORM REQUIREMENT R. 

C 

C INPUT DATA ARE THE MUMEER OF VERTICES, THE REQUIREMENT 
C R, AND THE NETWORK CONNECTION MATRIX. 

C 

C THE PROGRAM REAOS-IN THE INPUT DATA THEN DETERMINES THE 
C NUMBER OF VERTEX DISJOINT PATHS, STARTING WITH VERTEX 1 
C AS SOURCE, TO EACH OF THE REMAINING VERTICES. IT EMPLOYS 
C FRISCH'S ALGORITH.M FOR VERTEX-PAIR CONNECTIVITY AND 
C K LEI TRAM'S THEOREMS FOR OPTAINING THE OVERAIL CONNECT I- 
C V ITY. 

C 

C IF A CONNECTIVITY LESS THAN THE 'REQUIREMENT IS ENCOUNT- 
C ERED, THE PROGRAM CEASES FROM FURTHER VERIFICATION AND 
C PRINTS OUT THE CCNNECTIVITY VALUE SPECIFYING THE CORRESP- 
C. GNOING SOURCE VERTEX. OTHERWISE THE PROGRAM WILL INDICATE 
C THAT THE NETWORK IS AT LEAST R-CONNECTED. 

C 

C 

IMPLICIT INTEGER(A-Z} 

DIMENSION KMAT(50,S0J, ILISTUO) 

COMMON /ARR/ KMAT 
COMMON /CHAIN/ I LIST 

COMMON /PT/ SOURCE , SI NK,LL, G,N, RTEMP 
C 

READ (5,7J N,R 

READ (3,8) ( ( KMAT{ I , J) , J=1 ,30) t I ,N) 

L=0 

I NDEX=0 
RTEMP=R 
C 

DC 4 1=1, R 
IPLUS1=I+1 

DC 3 J=IPLUS1,N 



LIST DOWN PREVIOUSLY nESIGNAlED SINK VERTICES THAT 
ARE CONNECTED TO VERTEX J. 



LL = 0 

JLESSl=J-l 

IF ( JLESSl. LT. IPLUSl) GO TO 2 



1 



DC 1 L=IPLUS1, JL ESSl 
IF (KMAT{ L , J) .EQ.O) GO TO 1 



LL=LL+1 
ILIST (LL) = L 
CONTINUE 



IF (LL. GE.RTEMP) 60 TO 3 



2 



SEARCH FOR DISJOINT PATHS 
SOURCE=I 
SINK=J 



Q=0 

CALL VTXCON 
S=Q+LL+INDEX 
IF (S.GE.RJ 
GO TO 5 



GO TO 3 



ARE ALL (N-I) VERTICES PROCESSED? 
3 CONTINUE 



RTEMP=RTEMP-I 
INDEX = I NDEX-H 
4 CONTINUE 
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c 

c 

c 



c 



VERIFICATION OF RTH VERTEX AS SOURCE CCFPLETEU. 
WRITE (0,9) R 
GO TO 6 

CONNECTIVITY IS LESS THAN R. TERMINATE. 

5 WRl TE ( 6,10) S 
WRITE (6,11) I 

6 STOP 

7 FORMAT (216) 

8 FORMAT (2513/2513) 

9 FORMAT ( 'O' ,T10, ‘ NETWORK IS AT LE A ST • , I 3 , ' -CONNEC T EO ‘ ) 

10 FORMAT ('O', TIO, ‘NETWORK IS ONL Y *, 13, '-CONNEC TED WITH' 

11 FORMAT ('O', T13, * VE RT EX' , I 3 ,2 X , * AS SOURCE') 

END 



.... SUBROUTINE VTXCON ........... 

THIS SUBROUTINE DETERMINES THE f^UMbER OF VERTEX DIS- 
JOINT PATHS BETWEEN TWO GIVEN VERTICES USING FRISCH'S AL- 
GORITHM. IT IS HOWEVER MODIFIED BY I NCQRPORAT I KG KLEIT- 
(''.AN'S THEOREMS TO REDUCE COMPUTATIONAL TIME. 

SUBROUTINE VTXCON 
IMPLICIT INTEGER(A-Z) 

DlMfiNSlCN KMAT( 50,50), FLOMAT ( 50 , 50) 

DIMENSION NODE (50), J STORE (20), I STORE (20), I LI ST (10) 
CCMMON /ARR/ KMAT 
COMMON /CHAIN/ ILIST 

COMMON /PT/ SOUKCE,SINK,LL,Q,N, RTEMP 

INITIALIZE FLOW MATRIX TO ZERO 

DO 2 I=SOURCE,N 

DO 1 J=SOURCE,N 
FLOMAT ( I , J)=0 

1 CONTINUE 

2 CONTINUE 
MARK=0 

ERASE VERTEX LABELS 

3 DC A J=SOURCE,N 
N0DE(J)=0 

A CONTINUE 



INITIALIZE VARIABLES 
I = SOURCE 
L = 0 
M=0 
Z = 0 

NODE! I )=SOURCE 
SPLUS1=S0URCE+1 

IS SOURCE CONNECTED TO SINK? 

IF (KMAT (SOURCE, S INK) .EQ.O ) GO TO 5 
KMAT( SOURCE ,SINK) =0 
K MAT (SINK, SOURCE )=C 
Q = Q+1 
MARK=1 

SEARCH FOR AUGMENTING PATH 
5 DO 11 J=SPLUS1,N 

REMOVE PROCESSED VERTICES THAT ARE CONNECTED TO 
S INK 
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c 



IF (LL.EQ.O) GO TO 7 



6 



DO 6 K=1,LL 

I F (J.EO. ILISKKJ ) GO TO 11 
CONTINUE 



-7 ir CONNECTING PATH BETWEEN VERTEX-PAIR 

7 IF (FLOMAT ( J, I J . EO.O) GO TO 9 
IF (NODE! I) .LT.OJ GO TO 8 



IF 
I F 



A BACKWAfs.D BRANCH IS ENCOUNTERED 
(NODE! J J .LT.O) GO TO 11 
( J . bQ.SC.UPCE J GO TO 12 



8 

9 



VERTEX(J) IS EITHER UNLABELEC 
STRONGLY LABEL THIS VERTEX 
N0DE(J}=-I 

JSTORE(Mi=J 
GO TO 12 

IF (J. EQ. SOURCE) GO TO 11 



NODE! J)=-I 
GC TO 10 

I F (KMAT( I ,J) .NE.l) GO TO 11 



OR WEAKLY 



LABELEC. 



A FORWARD BRANCH IS ENCOUNTERED. TEST FOR ZERO FLOW 
IF (FLOMATI I » J) .NE.O) GO TO 11 
IF (N0DE( J ) .NE.O ) GO TO 11 
NODE! J) =I 



IS SINK VERTEX REACHED? 
IF (J.EU.SINK) GO TO 15 

10 M=f-H l 
JST0RE(M)=J 

11 CONTINUE 



12 IF (Z.NE.O) GO TO 14 
IF (M.EQ.O) GO TO 19 

TRANSFER CONTENTS OF JSTORE TO I STORE 

DO 13 S-1,M 
ISTORE(S)=JSTCRE(S) 

13 CONTINUE 



Z = M 
M = 0 

14 I-ISTORE(Z) 

Z=Z-1 

GO TO 5 

WE REACH HERE IF AN AUGMENTING PATH IS FOUND 
ADJUST ALL FLOWS ALONG THIS PATH. 

15 Q=Q+1 
U=G+LL 

IF (U.GE.RTEMP) GO TO 19 
TEST LABEL ON VERTEX J 

16 ]=N0DE(J) 

IF (I.LT.O) GO TO 17 

VERTEX IS WEAKLY LABELED. INCREMENT FORWARD FLOW BY 1 
FLOMATI I , J)=l 
GO TO 18 

VERTEX IS STRONGLY LABELED. DECREMENT BACKFLOW BY 1 

17 I=-I 
FLCMAT-( J, I J=0 

18 J = I 



IS SOURCE VERTEX REACHED? 
IF (J. EQ. SOURCE) GO TO 3 
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on 



GC TO 16 

VERTEX LABELING AND FLOW AUGMENTATION IS COMPLETED 

19 IF (MARK.EQ.O) GO TO 20 
KMAT( SOURCE , SINK) =1 
KHAT( SINK, SOURCE )=1 

20 RETURN 
END 
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